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ABSTRACT 


Selections were made for a conceptual plutonium 238 module for electric powers in 
the range from 2 to 10 kilowatts. Based on the selections, component and system effi- 
ciency and weight were estimated for powers of 2, 6, and 10 kilowatts. 
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ANALYSIS AND SELECTION OF DESIGN CONDITIONS FOR A 
RADIOISOTOPE BRAYTON-CYCLE SPACE POWERPLANT 
by John L. Klann 
Lewis Research Center 

SUMMARY 

A parametric analysis of a conceptual plutonium 238 radioisotope Brayton-cycle 
space powerplant was conducted. Major design variables were investigated to define 
design conditions which result in relatively compact and efficient power modules for 
electric powers in the range from 2 to 10 kilowatts. 

For the powerplant concept, a plane array of radioisotope capsules with a maximum 
capsule-wall temperature of 2460° R (1367° K) radiated its heat to the adjacent heat 
exchanger of the power conversion unit. In the power conversion unit, a mixture of 
helium and xenon drove a single combined rotating unit consisting of a single-stage, 
radial -flow compressor and turbine and a brushless generator of the Lundell type, all 
supported by gas bearings. Waste cycle heat was rejected by a separate liquid-filled 
coolant loop. 

As a result of this analysis, the major Brayton cycle parameter selections were: 
turbine inlet temperature, 2060° R (1144° K); compressor inlet temperature, 540° R 
(300° K); working fluid molecular weight, 83. 8; shaft speed, 36 000 rpm (3770 rad/sec); 
compressor pressure ratio, 1. 90; compressor specific speed, 0. 775; system loss pres- 
sure ratio, 0. 92; and recuperator effectiveness, 0.95. Based on these selections, over- 
all system efficiency (the ratio of net electrical power to initial isotope thermal power), 
unshielded system specific weights, and shielded specific weights were estimated as 
follows: 


Net 

power, 

kW 

Overall 

efficiency 

Unshielded 
specific weight 

Shielded 
specific weight 

lb/kWe 

kg/kWe 

lb/kWe 

kg/kWe 

2 

0. 18 

990 

449 

1590 

724 

6 

.23 

490 

222 

850 

386 

10 

.24 

440 

200 

740 

336 





INTRODUCTION 

Lewis Research Center is currently engaged in a program for advancing closed 
Brayton cycle technology with its ultimate intent of generating electric power in space. 
The initial efforts were centered about a solar -powered concept which is described in 
reference 1. The power conversion unit of that concept was designed to supply 8 kilo- 
watts of electric power at 400 hertz. It was characterized by two separate, main rotating 
units, namely, a turbine -driven compressor (turbocompressor) and a power -producing 
turboalternator. Some components for this concept are presently being tested. Although 
initially based on the solar -powered concept, the power conversion unit may also be 
employed with other energy sources. 

Reference 2 shows that radioisotopes are a promising heat source for a Brayton - 
cycle conversion unit. Successful development of safe, multikilowatt radioisotope heat 
sources coupled with Brayton- cycle power conversion modules would offer additional space- 
craft flexibility. Reference 3, for example, demonstrates the applicability of such a 
power system to a conceptual manned orbiting research laboratory. Two operative 5. 5- 
kilowatt (net electric) modules, each with an additional power conversion unit on standby, 
are used to enhance system reliability and maintenance. In comparison with solar- 
powered systems, the radioisotope -fueled systems do not require either orientation 
toward the Sun or large extended sunlight -collecting surfaces. 

The analysis presented herein was conducted to define desirable design conditions 
for an isotope -fueled power conversion module. Because of the limited availability of the 
more desirable radioisotopes (ref. 2), the basic goal was high overall system conversion 
efficiency. This goal was to be accomplished wherever possible with minimum weight 
but without excessive penalty to system volume. 

The analysis was restricted to a consideration of a single combined rotating unit - 
turbine, alternator, and compressor - supported by gas bearings and on a common shaft. 
The power range from 2 to 10 kilowatts was chosen for investigation. Although any given 
power system must have the capability of producing power on demand from zero to a 
nominal, rated value, the rated power output itself might also be adjusted for each given 
application. For a given Brayton combined rotating unit, this would be done by selecting 
a gas pressure level at, for instance, the compressor inlet. Hence, as part of this 
analysis, the feasibility of using one combined rotating unit over the range of rated power 
output from 2 to 10 kilowatts was also investigated. 

Gross shielding approximations and factors affecting fuel -block geometry and weight 
are presented in appendix D by John L. Klann, Leonard Soffer, and Gerald J. Barna. 

Brief descriptions of the calculations and typical results for each of the cycle heat- 
transfer components are presented in appendix E by John L. Klann, Gabriel N. Kaykaty, 
Paul T. Kerwin, and Dari D. Bien. 
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METHOD OF ANALYSIS 


In order to facilitate the analysis, a conceptual power system was hypothesized, and 
a reference set of conditions was established. Because of the large number of param- 
eters, each major variable was varied individually while the others were kept at the 
initial reference conditions. Thermodynamic variations were investigated first to deter- 
mine trends and desirable ranges of conditions for further analysis. Then, component - 
and system weight and volume variations were investigated. Based on this analysis of 
first-order variations, a set of preliminary operating conditions was selected. These 
conditions were then applied to the conceptual powerplant at 2-, 6-, and 10-kilowatt power 
levels, and overall system performance was evaluated. 


Conceptual Power System 

A schematic diagram of the assumed power conversion unit is shown in figure 1. 
Cycle state -point numbering is shown in figure 1. (All symbols are defined in appen- 
dix A.) The three major subsystems are a radioisotope heat source, a gas-filled power- 
conversion loop, and a liquid-filled heat-rejection loop. Radiation heat transfer between 
the heat source and the heat-source heat exchanger was assumed. The Brayton-cycle 
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Figure 1. - Schematic diagram of powerplant. 
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working fluid was assumed to be a mixture of helium and xenon (refs. 4 and 5). Recuper- 
ation was assumed. Both the compressor and turbine of the single-shaft combined rotat- 
ing unit were assumed to be single-stage radial -flow machines, while the alternator was 
a brushless generator of the Lundell type. The combined rotating unit was assumed to be 
supported hydrodynamically by bearings using the prime working gas. Cycle waste 
energy was assumed to be rejected by the secondary, liquid-filled coolant loop. This 
required a gas-to-liquid waste heat exchanger, a radiator, and a pump. The pump was 
assumed to be motor driven and separate from the prime turbomachinery. 


Thermodynamic Calculations 

Single-stage radial -flow turbine and compressor efficiency models were hypothesized 
and empirically fitted to experimental results obtained from component testing (ref. 6). 
Efficiency was assumed to be a function of specific speed, specific diameter, Reynolds 
number, characteristic tip Mach number, and tip diameter. These models (presented in 
appendix B) were included in the basic cycle thermodynamic calculation procedure 
(presented in appendix C). Hence, effects of machinery size, speed, and pressure level 
on cycle thermal efficiency (the ratio of gross shaft power to cycle thermal input power 
Psh/Qin) were evaluated. 

Table I shows the initial reference conditions and the range of parameters investi- 
gated for basic cycle performance. The helium-xenon gas mixture (used because of its 
increased heat -transfer capability in comparison with pure gases) was investigated for 
the average molecular weights of krypton (83. 8) and also, but less thoroughly, of argon 
(39. 9). Although the use of a gas mixture is advantageous, the ability to maintain the 
proper molecular weight remains to be demonstrated. Limiting the analysis and design 
considerations to the molecular weights of krypton and argon, rather than to a range of 
molecular weights, would permit substitution of the single gas without an appreciable 
effect on rotating component efficiency. Such a substitution, however, would result in 
reduced overall conversion efficiency because of reduced heat transfer. 

Turbine inlet temperatures investigated were 1860°, 1960°, and 2060° R (1033°, 

1089°, and 1144° K). A basic guideline for the analysis was to use nonrefractory metals 
for all components except the heat source. The selected range of turbine inlet tempera- 
ture was presumed to be permissible for long-term use of nonrefractory alloys. 

At each turbine inlet temperature, compressor inlet temperature was varied over a 
range bounded by 490° and 720° R (272° and 400° K). The other major variables included 
recuperator effectiveness, system loss pressure ratio (i. e. , that portion of the compres- 
sor pressure ratio which is available to the turbine), compressor pressure ratio, shaft 
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TABLE I. - INITIAL REFERENCE CONDITIONS AND RANGE OF 


PARAMETERS FOR THERMODYNAMIC CALCULATIONS 


Parameter 

Initial 

reference 

condition 

Range or 
other 
conditions 

Working fluid 

Helium -xenon 


Molecular weight, M w 

83.8 

39.94 

Turbine inlet temperature, T.«: 
°R 

1960 

1860, 2060 

°K 

1089 

1033, 1144 

Compressor inlet temperature, T 4 : 
°R 

540 

490 to 720 

°K 

300 

272 to 400 

Recuperator effectiveness, E r 

0.900 

0. 875 to 0.950 

Loss pressure ratio, L 

0.900 

0.900 to 0.950 

Compressor pressure ratio, P 5 /P 4 

1.90 

1.60 to 3.60 

Rotational speed, N: 
rpm 

36 000 

24 000; 48 000; 60 000 

rad/ sec 

3770 

2510; 5020; 6280 

Compressor specific speed, N g c 

0.775 

0. 54 to 0. 86 

Turbine -exit -loss parameter, L g 

0. 04 


Equivalent sink temperature, T g : 
°R 

450 


°K 

250 


Radiator surface emittance, e 

0.90 


Pres sure -drop distribution ratio, ot 

0. 50 


Compressor slip factor, j3 

0. 85 


Turbine blade -jet speed ratio, v 

0.70 


Waste -heat -exchanger effectiveness, E g 

0.95 


Waste -heat -exchanger capacity rate ratio, C g 

0. 87 



rotational speed, and compressor specific speed. The rotational speeds investigated 
were chosen to yield multiples of 400-hertz output frequencies for four -pole alternators. 

The remaining parameters (listed in table I) were assigned reasonable values but 
were not varied. The turbine -exit -loss parameter, defined as the square of the ratio of 
the axial component of the velocity to the ideal jet speed, was set at 0. 04. Reference 7 
shows this to be a lower limit for stable flow conditions in radial inflow turbines. This 
loss parameter is a measure of the turbine -exit kinetic energy and relates total to static 
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efficiency (eq. (C24)). Turbine static efficiency has been used in this analysis, and hence 
the exit kinetic energy, although made small, was assumed to be lost to the system. The 
assumed equivalent sink temperature is a few degrees below the maximum equivalent 
sink temperature for an unoriented, cylindrical radiator in a low Earth orbit. The 
assumed radiator surface emittance is within the state of the art for low -temperature 
surfaces. The pressure -drop distribution ratio is the ratio of fractional pressure drop 
Ap/p on the high-pressure side of the cycle to that on the low-pressure side. The com- 
pressor slip factor and the turbine blade to jet speed were taken as representative values 
for efficient designs. 


Weight Calculations 

Results of the thermodynamic calculations were used as input conditions for the 
determination of size and weight tradeoffs among the components. The same approach 
used in the thermodynamic calculations was used for the weight calculations. Each major 
variable was investigated individually while the others were kept at the initial reference 
conditions. The variables were turbine inlet temperature, compressor inlet temper- 
ature, recuperator effectiveness, system loss pressure ratio, and to a limited extent, 
working fluid and its molecular weight. Shaft power level was set at 6. 875 kilowatts 
as representative of the midrange of net system output. 

In addition, particular geometries were assumed for each component. A manned 
5-year mission was postulated for sizing the heat source and radiator. As in reference 2, 
the heat source was assumed to be a flat slab (referred to as "fuel block") formed by a 
matrix of closely packed fuel capsules joined together by a niobium web. The fuel was 
assumed to be plutonium 238 in a dioxide form. The fuel block was sized to be capable 
of rejecting its initial heat production without normally exceeding a temperature of 
2460° R (1367° K). The initial inventory allowed for isotope decay over the assumed 
5-year mission and an assumed 90-day launch delay. Single -sided thermal radiation from 
one face of the fuel block to the heat-source heat exchanger was assumed. 

Figure 2 depicts the assumed installation of the power system in a mission vehicle. 
The power system was assumed to be placed near the cylindrical wall at the aft end of a 
260-inch- (6. 6-m-) diameter manned vehicle. The radiator was assumed to encompass 
this portion of the vehicle. Immediately inside the thermally insulated door is the radio- 
isotope heat source mounted on a reentry body. Nuclear shielding was assumed to 
enclose the fuel block and heat-source heat exchanger on all sides excepting those which 
face away from the manned compartment. Thermal insulation might be required on the 
inside of the shield in order to keep the shielding material at reasonable temperatures. 
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The hot gas mixture is ducted from the heat-source heat exchanger around the shield to 
the power conversion unit. 

Analysis of the fuel block and gross approximation of the shield thickness and weight 
are presented in appendix D. For the general analysis of fuel block size and weight, one 
layer of fuel capsules was assumed. Each capsule had a void volume to allow for a maxi- 

n 

mum internal helium pressure of 10 000 psi (69 MN/m ) (see ref. 2 for details). Capsule 
burial criteria were not considered. Shield weights were calculated for constant nuclear 
dose rates of 1 and 5 millirems per hour at a distance of 10 feet (3 m) from the source. 
These dose rates were considered to be representative of conservative and liberal allow- 
ances, respectively, for a manned spacecraft. 

Details of the four heat-exchanger analyses are presented in appendix E. The heat- 
source heat exchanger was assumed to have a once-through plate-fin arrangement, while 
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the recuperator employed a rectangular off set -fin counterflow core with triangular end 
sections. The heat-sink heat exchanger was assumed to have a multipass cross counter- 
flow tube -fin arrangement. This exchanger was sized to fit directly on the appropriate 
recuperator end section. A cylindrical, unoriented radiator with longitudinal finned tubes 
and circumferential headers on the 260-inch- (6. 6-m) diameter vehicle was assumed. 

Two complete coolant circuits (one redundant) sharing common fins were so armored 
against meteroid penetration that the probability of survival of at least one coolant circuit 
at the end of 5 years was 0. 995. The fluorocarbon FC-75 (ref. 8, p. 86) was selected 
as the radiator coolant. 


System Efficiency 

In the evaluation of overall system efficiency, calculations or allowances were made 
for some of the thermal, mechanical, and electrical losses. Details are presented in 
appendix F. Allowances were made for heat-source thermal loss, system control power, 
and liquid -coolant -loop pump power. Estimates were made for bearing friction loss, 
alternator electromagnetic and windage losses, and shaft seal loss for the combined 
rotating unit. No direct allowance was made for thermal loss through the insulation 
surrounding the conversion loop. Although it was assumed that the combined rotating 
unit alternator and the control system electronics would be cooled by the radiator coolant, 
these heat loads were not included in the radiator calculations. These heat loads might 
add 15 to 25 percent to the rejected heat for power outputs of 10 to 2 kilowatts, respec- 
tively, but because these heat loads do not change rapidly with the design variables in 
table I, their inclusion in or omission from the analysis has only a minor effect on selec- 
tion of desirable design conditions for the power system . 

’’System efficiency” is herein defined as the ratio of alternator output power less 
pump and control power to the gross initial thermal energy from the isotope package. 
Gross initial isotope energy allowed for decay over the mission life, heat-source thermal 
loss, and required cycle thermal input, and no credit for the initial excess thermal 
energy available from the heat source was taken. 


System Variations with Power Level 

Design parameters from this first-order thermodynamic and weight analysis were 
used to define the conceptual powerplant. Its efficiency and total weight were estimated 
for net output power levels of 2, 6, and 10 kilowatts. The combined rotating unit was 
sized for operation over the entire power range, while all other components were resized 
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at each power level. Hence, the selected values for turbine inlet temperature, compres- 
sor inlet temperature, recuperator effectiveness, and system loss pressure ratio were 
maintained at each power level. Although the system loss pressure ratio was fixed, the 
individual conversion-loop-component pressure drops were redistributed at each power 
level in order to yield minimum combined heat -exchanger weight (appendix E). 

Total system weight estimates are based on the power system arrangement in 
figure 2. The heat-source subsystem weight estimates include the fuel block, reentry 
body, and a quantity of emergency heat-storage material. This combined heat-source 
package is configured to contain the radioisotope fuel safely under a variety of conditions. 
Details of the concept are presented in reference 2 and further explored in references 9 
and 10. The reentry body permits separate intact return of the isotope capsules from 
orbit, while the main purpose of the heat-storage material is to prevent capsule melting 
and isotope release as the result of a launch-pad explosion and fire and subsequent entrap- 
ment of the heat source beneath the burned out wreckage of the launch vehicle. In the present 
analysis, enough heat- storage material was assumed (based on the work of ref. 10)for approx- 
imately a 90-minute delay in recovery of the heat source after such an accident. 

On the basis of minimum combined fuel -block and shield weight, the fuel block was 
assumed to have one layer of capsules with each capsule void volume sized for a maxi- 
mum internal pressure of 12 500 psi (86. 1 MN/m^) rather than the 10 000 psi (69 MN/m^) 
assumed in Weight Calculations which resulted in somewhat higher weight (appendix D). 
Both the reentry body and heat-storage weights were estimated from the results of refer- 
ences 9 and 10. The reentry body was assumed to be approximately 20 percent of the 
total heat-source weight, while the heat-storage material was assumed to be beryllium 
oxide estimated from references 9 and 10 as 150 pounds (68 kg) for a 25-kilowatt-thermal 
source. 

The heat -exchanger weights were estimated from the analysis in appendix E. The 
weight of the combined rotating unit was based on a detailed design. Ducting, structure, 
and control system weights were scaled from the results of reference 3. A specific 
radiator area of 90 square feet per kilowatt electric (8. 36 mVkWe) was assumed so that 
radiator weight near the minimum could be achieved at each power level. Structural 
aspects of the radiators were not considered. The shield weight estimates were based on 
a dose rate of 5 millirems per hour at a distance of 10 feet (3 m) from the isotope heat 
source. 


RESULTS AND DISCUSSION 

First-order variations of major power -conversion -unit parameters are presented 
in the sections Thermodynamic Variations and Weight Variations. The final section 
Selections and Performance Estimates summarizes the results. 
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Thermodynamic Variations 


A performance map for the initial reference conditions of table I is shown in figure 3. 
Cycle thermal efficiency (the ratio of gross shaft power to required cycle thermal input, 
where gross shaft power is the turbine power minus compressor power) is plotted 
against specific prime radiator area (the ratio of prime or minimum radiator area 1 to 
gross shaft power) with compressor pressure ratio and inlet temperature as independent 
parameters. For a fixed compressor pressure ratio, cycle efficiency approaches a 
maximum as the compressor inlet temperature approaches the sink temperature (450° R, 
or 250° K). At this condition, the largest cycle temperature difference, and hence 
efficiency, is available; however, the radiator area becomes large. For a fixed value of 
compressor inlet temperature, cycle efficiency exhibits a maximum with respect to 



4 6 8 10 12 

Specific prime radiator area, A r p / P sh , irr/kW 


Figure 3. - Brayton cycle performance. (All other initial reference 
conditions of table I apply.) 


1 In the thermodynamic calculations, prime or minimum radiator area is approxi- 
mated by neglecting the drop between the bulk liquid temperature and the radiating sur- 
face temperature (eq. (C36)). 
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compressor pressure ratio. The variations in specific prime radiator area with com- 
pressor inlet temperature and pressure ratio result from combined effects of radiator 
heat load and average radiator temperature level. 

The general level of performance shown in figure 3 is dependent not only on the 
assigned cycle conditions, but also on the calculated turbine and compressor efficiencies 
(appendix B). The ranges of efficiency calculated for this figure were 0. 847 to 0. 855 
and 0. 754 to 0. 771 for the turbines and compressors, respectively. Since these varia- 
tions were relatively small, the general shape of these curves is similar to that reported 
in other analyses using fixed turbine and compressor efficiencies (ref. 11, e.g. ). 

Turbine- and compressor -inlet temperatures. - Figure 4 shows the variation of 
cycle thermal efficiency and specific prime radiator area with compressor inlet temper- 
ature for three turbine inlet temperatures. Compressor inlet temperature is shown 
only for the range of interest, namely, 500° to 630° R (278° to 350° K), where the high 
value is close to that for minimum specific prime radiator area. At the reference value 
of compressor inlet temperature (540° R, or 300° K), cycle efficiency increased by 
approximately 0. 02 (or about 8 percent) for each 100° R (56° K) increase in turbine 




280 300 320 340 360 

Compressor inlet temperature, °K 

Figure 4. - Effect of compressor- and turbine-inlet 
temperature on specific prime radiator area and 
cycle thermal efficiency. (Ail other initial 
reference conditions of table I apply.) 
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inlet temperature above 1860° R (1033° K), but specific prime radiator area decreased 
by approximately 6. 5 square foot per kilowatt or 0. 6 square meter per kilowatt (or about 
10 percent) for the same changes. For this example, predicted turbine and compressor 
efficiencies decreased from 0. 850 and 0. 765 to 0. 847 and 0. 763, respectively, between 
1860° and 2060° R (1033° and 1144° K). 

On the basis of cycle performance alone, the highest tolerable turbine inlet temper- 
ature should be chosen. 

Compressor specific speed. - Variations in turbine- and compressor -efficiency, 
size, and pressure level due to changes in compressor specific speed are a direct result 
of the assumed turbomachinery efficiency models (see appendix B). Variation in cycle 
efficiency and specific prime radiator area with compressor pressure ratio and specific 
speed is shown in figure 5(a). As compressor specific speed increases from 0. 5425 to 
0.7750 for any value of compressor pressure ratio, cycle efficiency increases, and 
specific prime radiator area decreases. Near -maximum cycle efficiencies are obtained 
for compressor specific speeds from 0. 6975 to 0. 7750. Increasing compressor specific 
speed from 0.7750 to 0. 8525 reduces cycle efficiency. Maximum cycle efficiency occurs 
at a compressor pressure ratio of about 2. 1 for any compressor specific speed. 

Maximum system efficiency tends to occur at slightly lower values of compressor 
pressure ratio than those which result in maximum cycle efficiency. Among the parasitic 
losses included in the calculation of system efficiency, alternator windage decreases and 
bearing and seal losses increase with increasing compressor pressure ratio, whereas 
the others remain essentially constant (see appendix F for details). 

Figure 5(b) shows the effect of compressor specific speed on turbine and compressor 
efficiencies, the specific compressor inlet pressure (the ratio of compressor inlet 
pressure to gross shaft power), and the required tip diameters, all over a range of com- 
pressor pressure ratio from 1. 5 to 2. 8. Decreasing compressor pressure ratio in- 
creases the specific compressor inlet pressure and decreases turbine and compressor 
diameters. At the same time, compressor efficiency improves (mainly because of the 
beneficial effect of pressure ratio; see appendix B) and turbine efficiency remains 
essentially unchanged (because of conflicting effects of size, pressure ratio, and 
Reynolds number; see appendix B). For any compressor pressure ratio, decreasing 
compressor specific speed increases the specific compressor inlet pressure, slightly 
increases compressor tip diameter, and decreases compressor efficiency. Turbine 
efficiency, however, first increases and then decreases with decreasing compressor 
specific speed, when the compressor pressure ratio is fixed. 

One reason for selecting a compressor specific speed below 0. 775 might be to 
increase overall system pressure level. For example, at the reference pressure ratio 
(1.90), reductions of compressor specific speed from 0.7750 to 0.6975, 0.6200, and 
0. 5425 increased compressor inlet pressure (for a fixed shaft power) by about 23, 62, 
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4 5 6 7 8 9 

Specific prime radiator area, A r p / P sh , mWkW 

(a) Specific prime radiator area and cycle thermal efficiency 
variations. 




Compressor 
specific speed, 
N s,c 

a-0.6200, .6975 


.5425, .7750 


: } Turbine 



Compressor 



5 1 — 



Compressor pressure ratio 


(b) Pressure level and turbomachinery variations. 


Figure 5. - Effect of compressor specific speed on performance. (All other initial reference conditions of table I apply. ) 
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Specific prime radiator area, A rp /P shp fr/kW 


4 5 6 7 8 2 9 

Specific prime radiator area, A rp /P sh , nr/kW 

(a) Specific prime radiator area and cycle thermal efficiency 
variations. 
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(b) Pressure level and turbomachinery variations. 


Figure 6. - Effect of turbomachinery rotational speed on performance. (All other initial reference conditions of table 1 apply). 
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and 150 percent, respectively. Since one design goal for this analysis was to consider 
use of one combined rotating unit over the net system power range from 2 to 10 kilowatts, 
it was required that both compressor specific speed and pressure ratio be selected so 
that reasonably high bearing -cavity pressures were available to the gas bearings of the 
rotating unit at the 2 -kilowatt power level. 

Rotational speed. - Results of the basic thermodynamic cycle calculations are 
presented in figure 6 for three rotational speeds. These speeds were chosen to be com- 
patible with the supposition (see METHOD OF ANALYSIS) that only multiples of 
400 hertz should be considered. Hence, the rotational speeds 24 000, 36 000, and 
48 000 rpm (2510, 3770, and 5020 rad/sec) coupled to a four-pole alternator produce, 
respectively, 800, 1200, and 1600 hertz. The performance variations due to rotational 
speed are a direct result of the assumed turbomachinery efficiency models. 

Figure 6(a) shows that cycle performance improves as rotational speed is reduced. 
Peak cycle efficiencies (at a compressor pressure ratio of about 2. 1) show a 0. Oil 
(4 percent) increase as rotational speed is reduced from 48 000 to 36 000 rpm (5020 to 
3770 rad/sec). Reducing rotational speed from 36 000 to 24 000 rpm (3770 to 2510 
rad/sec) increases cycle efficiency by 0. 013 (4. 6 percent). Figure 6(a) presents the 
effect of rotational speed, over a range of compressor pressure ratio, on turbine and 
compressor efficiencies and tip diameters as well as specific compressor inlet pressure. 
Decreasing rotational speed necessitates increasing tip diameters but results in increas- 
ing compressor efficiency, while system pressure levels decrease. Reducing shaft 
speed from 48 000 to 36 000 rpm (5020 to 3770 rad/sec) improves turbine efficiency only 
slightly over the range of compressor pressure ratio investigated, and the change from 
36 000 to 24 000 rpm (3770 to 2510 rad/sec) improves turbine efficiency only at low 
pressure ratios (below 2.1). At the lower shaft speed, the reduced turbine efficiency 
improvement results from the size correlation employed in appendix B. In appendix B, 
it was assumed that size effects occur for radial turbines and compressors below a 
6-inch (15-cm) tip diameter. As may be observed, the turbines and compressors at 
24 000 rpm (2510 rad/sec) exceed 6-inch (15-cm) tip diameters above compressors 
pressure ratios of 1. 62 and 2. 03, respectively. 

Choice among these speeds would be made on the basis of cycle efficiency, machinery 
size, and system pressure level. Although the 24 000-rpm (2510-rad/sec) solutions 
exhibit higher cycle efficiencies, low system pressure level might be a problem. Of 
course, the possibility of reducing compressor specific speed below 0. 7750 tends to 
alleviate this drawback. However, it should also be noted that lower compressor pres- 
sure ratios are required to attain higher specific compressor inlet pressures. There is 
a tendency for a Brayton-cycle system to become more sensitive to temperature and/or 
pressure variations at low values of machinery pressure ratios. Decreasing compressor 
pressure ratio decreases compressor and turbine work and the temperature differences 
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across these components. Hence, at low compressor pressure ratio, small changes in 
the inlet temperatures result in a larger change in the component work and might more 
greatly affect the system dynamics. The solutions for the 48 000-rpm (5020-rad/sec) 
speed exhibit higher pressure levels but lower cycle efficiencies and smaller turbine 
and compressor tip diameters. 

Recuperator effectiv eness. - Typical effects of recuperator effectiveness on cycle 
performance are shown in figure 7. A range of recuperator effectiveness and compres- 
sor pressure ratio are plotted to show their overall effects on cycle thermal efficiency 
and specific prime radiator area. In general, increasing recuperator effectiveness 
increases cycle efficiency and decreases specific prime radiator area. Although in- 
creasing recuperator effectiveness decreases the average radiator temperature, the 
increased cycle efficiency reduces the radiator heat load, and the net result is a decrease 
in required radiator area. 

The changes in maximum cycle efficiency with changes in recuperator effectiveness 
from 0. 875 to 0. 900, 0. 900 to 0. 925, and 0. 925 to 0. 950 were 0. 011 (4. 0 percent), 

0. 13 (4.6 percent), and 0.017 (5.7 percent), respectively. Hence, for equal increments 
in effectiveness, cycle efficiency increases more rapidly as effectiveness itself rises. 

It should be noted that the compressor pressure ratio required to produce maximum 
cycle efficiency at each recuperator effectiveness decreases from approximately 2. 2 to 



4 5 6 7 8 9 

Specific prime radiator area, A r p/P^, nr/ kW 

Figure 7. - Effect of recuperator effectiveness on performance. 
(All other initial reference conditions of table I apply.) 
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1. 9 as recuperator effectiveness is increased from 0. 875 to 0. 950. In this respect, the 
level of recuperation will affect the selection of compressor pressure ratio and hence the 
machinery size and system pressure levels. 

System loss pressure ratio . - Figure 8 shows the effects on cycle thermal efficiency 
and specific prime radiator area from varying system loss pressure ratio and compres- 
sor pressure ratio. At any value of compressor pressure ratio, increasing the system 
loss pressure ratio (decreasing system pressure losses) increases cycle efficiency and 
decreases specific radiator area. 

The changes in peak system efficiency from changes in loss pressure ratio from 
0. 900 to 0. 925 and from 0. 925 to 0. 950 were 0. 023 (8. 1 percent) and 0. 021 (6. 9 percent), 
respectively. Hence, for equal -increment increases in loss pressure ratio, the increases 
in cycle efficiency decrease as system losses are reduced. Those compressor pressure 
ratios that produce maximum cycle efficiency at each loss pressure ratio decrease from 

2. 1 to 1. 9 with increasing system loss pressure ratio from 0. 900 to 0. 950. At any value 
of compressor pressure ratio, increasing loss pressure ratio has a negligible effect on 
turbine and compressor efficiencies in spite of decreases in wheel diameters and system 
pressure level. For example, at a compressor pressure ratio of 1.9, turbine and com- 
pressor efficiencies decrease by 0. 001 and 0. 002 in going from a loss pressure ratio of 
0. 900 to 0. 950, respectively, while the specific compressor inlet pressure decreases 


Compressor <8> Initial reference 

pressure condition 



24 40 56 72 88 

Specific prime radiator area, A r p /P sh , ft 2 /kW 


3 4 5 6 7 8 

Specific prime radiator area, A r p/P^, m 2 /kW 


Figure 8. - Effect of loss pressure ratio on performance. (All 
other initial reference conditions of table I apply.) 
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Cycle thermal efficiency, p sh ,Q ii 
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4 5 6 7 8 

Specific prime radiator area, A r p/p sh , nr/kW 

(a) Effect of rotational speed. (b) Pressure level and turbomachinery variations. 

Figure 9. - Performance at working fluid molecular weight of 39. 9. (All other initial reference conditions of table I apply. ) 
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from 2. 0 to 1.6 psia per kilowatt or 13. 8 to 11 newtons per meter squared kilowatt (a 
20-percent decrease). 

Work ing -fluid molecular weight . - Effects on cycle performance of the helium - 
xenon mixture at the molecular weight of argon are explored in figure 9 (similar to fig. 6 
for molecular weight of krypton) for three rotational speeds. Figure 9(a) shows that 
peak cycle efficiencies occur at a compressor pressure ratio of about 2. 1. Decreasing 
rotational speed from 60 000 to 48 000 rpm (6280 to 5020 rad/sec) increases cycle effi- 
ciency 0. 008 (3 percent), and a decrease in speed from 48 000 to 36 000 rpm (5020 to 
3770 rad/sec) increases cycle efficiency 0. 010 (3. 7 percent). Hence, at this lower 
molecular weight, the sensitivity of cycle efficiency to rotational speed was about the 
same as at the higher molecular weight. 

Figure 9(b) shows the effect of changing rotational speed on turbine- and compressor- 
tip diameters, component efficiencies, and specific compressor inlet pressure over a 
range of compressor pressure ratio. The trends are similar to those for rotational 
speeds of 24 000, 36 000, and 48 000 rpm (2510, 3770, and 5020 rad/sec) at the molecu- 
lar weight of 83. 8. Because the specific compressor inlet pressure must be chosen high 
enough to permit successful bearing operation at the 2-kilowatt-electric system power 
level, specific compressor inlet pressures are greatly reduced from those at the same 
speed but higher molecular weight, an effect that is important. For example, for 36 000- 
rpm (3770-rad/sec) solutions and a compressor pressure ratio of 1.9, compressor spec- 
ific inlet pressure is reduced about 65 percent from approximately 2. 08 psia per kilowatt 
(14. 3 (kN/m^)/kW) at the 83. 8 molecular weight to 0. 72 psia per kilowatt (4. 96 
(kN/m /kW) at the 39. 9 molecular weight. Furthermore, specific compressor inlet pres- 
sures for the 36 000-rpm (3770-rad/sec) solutions at 83. 8 molecular weight are nearly 
the same as those for 60 000-rpm (6280-rad/sec) solutions at 39.9 molecular weight. 

Figure 10 compares the results for 36 000 rpm (3770 rad/ sec) and the molecular weight 
of 83. 8 with the results for 60 000 rpm (6280 rad/sec) and the molecular weight of 39. 9. 
Over the range of compressor pressure ratio investigated, the higher molecular weight 
solutions have higher cycle efficiencies, lower prime radiator areas, and larger turbine 
and compressor tip diameters. In addition, maximum cycle efficiency shows an increase 
of 0. 020 and specific radiator area shows a decrease of 6. 25 square feet per kilowatt 

Q 

(0. 58 m /kW) for the higher molecular weight. In addition to these apparent performance 
advantages, it was also calculated that, for the higher molecular weight, the turbine-exit 
impeller -root stress would be approximately half that for the lower molecular weight; this 
effect could be important in trying to meet the design goal of a 5-year minimum life for 
the system. 
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TABLE H. - CYCLE CONDITIONS FOR SIZE AND WEIGHT COMPARISONS 


[Shaft power, 6 . 875 kW. ] 


(a) U.S. Customary Units 


Parameter a 

Turbine inlet 
temperature, 

> 

°R 

Compressor inlet 
temperature, 

T 4> 

°R 

Recuperator effectiveness, 
E 

'Loss pressure ratio, 
L 

Molecular weight, 
M w 

1860 

I960 

(b) 

2060 

490 

540 

(b) 

588 

0. 875 

0.900 
0 b) 

0. 925 

0.950 

0.900 

(b) 

0. 925 

0.950 

C 39. 9 

83.8 

(b) 

Turbine efficiency, 

0.850 

0. 849 

0.847 

0.847 

0.849 

0.850 

0.849 

0.849 

0.849 

0. 849 

0. 849 

0. 849 

0,848 

0.835 

0.849 

Compressor efficiency, 7] c 

0.765 

0.764 

0.763 

0.762 

0.764 

0.759 

0.764 

0.764 

0.764 

0.764 

0.764 

0.763 

0.762 

0.751 

0.764 

Cycle efficiency, 77 

0.257 

0.279 

0.299 

0.317 

0.279 

0.235 

0.265 

0.279 

0.295 

0.313 

0.279 

0.304 

0.327 

0.261 

0.279 

2 

Prime area, A^ . ft 
r, p> 

423 

372 

332 

628 

372 

292 

385 

372 

358 

343 

372 

330 

297 

404 

372 

Mass flow rate, m, lb/ sec 

1. 109 

0.960 

0. 847 

0.832 

0.960 

1 . 160 

0.960 

0.960 

0.960 

0. 960 

0. 960 

0.853 

0.770 

0. 495 

0. 960 

Temperature, °R: 
















Turbine exit, T 2 

1555 

1639 

1723 

1639 

1639 

1638 

1639 

1639 

1639 

1639 

1639 

1624 

1610 

1644 

1639 

Heat-sink heat -exchanger inlet, T 3 

826 

835 

844 

775 

835 

897 

857 

835 

813 

790 

835 

834 

833 

839 

835 

Compressor exit, Tg 

745 

746 

746 

678 

746 

815 

746 

746 

746 

746 

746 

746 

746 

749 

746 

Heat-source heat -exchanger inlet, T g 

1474 

1549 

1625 

1543 

1549 

1556 

1527 

1549 

1572 

1594 

1549 

1536 

1524 

1554 

1549 

Pressure, psia: 
















Turbine inlet, p 1 

29.4 

25.4 

22.4 

23.1 

25.4 

29.4 

25.4 

25.4 

25.4 

25.4 

25.4 

22.8 

20.8 

25.1 

25.4 

Turbine exit, p 2 

17.2 

14.9 

13.1 

13.5 

14.9 

17.2 

14.9 

14.9 

14.9 

14.9 

14.9 

13.0 

11.5 

14.7 

14.9 

Compressor inlet, p^ 

16.0 

13.8 

12.2 

12.6 

13.8 

16.0 

13.8 

13.8 

13.8 

13.8 

13.8 

12.3 

11 . 1 

13.7 

13.8 

Compressor exit, p^ 

30.4 

26.3 

23.2 

23. 9 | 

26. 3 j 

30.4 

26.3 

26.3 

26.3 

26.3 

26.3 

23.4 

21.1 

26.0 

26.3 


a At reference conditions of table I unless otherwise noted. 


^Initial reference conditions. 
c Rotational speed, 60 000 rpm. 


TABLE H. - Concluded. CYCLE CONDITIONS FOR SIZE AND WEIGHT COMPARISONS 


[Shalt power, 6. 875 kW. ] 


(b) SI Units 


Parameter 3- 

Turbine inlet 
temperature, 

> 

°K 

Compressor inlet 
temperature, 

> 

°K 

Recuperator effectiveness, 
E 

Loss pressure ratio, 
L 

Molecular weight, 
M w 


j 1033 

1089 

(b) 

1144 

; 272 

300 

(b) 

327 

0. 875 

0.900 

(b) 

0.925 

0.950 

0.900 

(b) 

0.925 

0.950 

c 39.9 

83.8 

(b) 

Turbine efficiency, 

0. 850 

0.849 

0.847 

0. 847 

0.849 

0. 850 

0.849 

0.849 

0. 849 

0. 849 

0. 849 

0. 849 

0.848 

0.835 

0.849 

Compressor efficiency, ?7 C 

0.765 

0.764 

0.763 

0.762 

0.764 

0.759 

0.764 

0.764 

0.764 

0.764 

0.764 

0.763 

0.762 

0.751 

0.764 

Cycle efficiency, r j 

, 0.257 

0.279 

0.299 

0.317 

0.279 

0.235 

0.265 

0.279 

0,295 

0.313 

0.279 

0.304 

0.327 

0.261 

0.279 

„ . . 2 
Prime area, A f p , m 

40.2 

34.6 

30.8 

58.3 

34.6 

27. 1 

35.8 

34.6 

33.3 

31.8 

34.6 

30.6 

27.6 

37.5 

34.6 

Mass flow rate, m, kg/ sec 

0.459 

0.436 

0.385 

0. 378 

0. 436 

0. 526 

0. 436 

0. 436 

0. 436 

0.436 

0. 436 

0. 387 

0.350 

0.225 

0. 436 

Temperature, °K: 
















Turbine exit, T 2 

864 

910 

957 

910 

910 

909 

910 

910 

910 

910 

910 

903 

894 

914 

910 

Heat- sink heat- exchanger inlet, Tg 

459 

464 

468 

430 

464 

498 

476 

464 

452 

439 

464 

463 

462 

466 

464 

Compressor exit, Tg 

419 

420 

420 

376 

420 

452 

420 

420 

420 

420 

420 

420 

420 

416 

420 

Heat -source heat-exchanger inlet, Tg 

819 

860 

902 

858 

860 

864 

848 

860 

874 

886 

860 

853 

847 

864 

860 

Pressure, kN/m : 

Turbine inlet, p 1 

203 

175 

154 

159 

175 

203 

175 

175 

175 

175 

175 

157 

143 

173 

175 

Turbine exit, pg 

119 

103 

90 

93 

103 

119 

103 

103 

103 

103 

103 

90 

79 

101 

103 

Compressor inlet, p^ 

110 

95 

84 

87 

95 

110 

95 

95 

95 

95 

95 

85 

77 

94 

95 

Compressor exit, pg 

210 

181 

160 

165 

181 

210 

181 

181 

181 

181 

181 

161 

145 

179 

181 


a At reference conditions of table I unless otherwise noted. 


^Initial reference conditions. 
c Rotational speed, 6280 rad/ sec. 



Weight Variations 


Table II shows the typical cycle conditions which were used for the calculations of 
component and system weights. For each parameter varied, as shown in table n, the 
remaining cycle parameters were held at the initial reference values of table I. The 
shaft power was assumed to be 6. 875 kilowatts. 

It is recognized that the results of this section are not completely general. Particu- 
lar mission requirements or other component models might alter these results. For 
example, more compact core geometries could reduce heat -exchanger weights, life- 



Figure 11. - Effect of turbine inlet temperature on major 
component weights. Gross shaft power, 6.875 kilowatts, 
(All other initial reference conditions of table I apply. ) 
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support -system thermal requirements could reduce the power -system -radiator heat 
loads, existing vehicle structure could be used for the radiator and its weight not be 
charged to the power system, heat-source geometry and shielding requirements could be 
varied considerably, etc. Of all these possibilities, none appears to have considerable 
influence on this first-order selection of Brayton-cycle parameters. The assumed 
component models are considered to be representative of current thinking and form a 
reasonable basis for the preliminary selections of turbine inlet temperature, compressor 
inlet temperature, recuperator effectiveness, and system loss pressure ratio. 

Turbine inlet temperature. - Figure 11 presents major component weights and their 

sum as a function of turbine inlet temperature. Radiator, fuel -block, and shield weights 

are presented. The radiator area has been assumed to be 650 square meters (approxi- 
2 2 

mately 95 ft /kW or 8. 8 m /kW). This area resulted in near -minimum weight for the 
radiator over the range of turbine inlet temperature. Bands of shield weights and corre- 
sponding total weights are shown bounded by lines of constant nuclear dose rate of 1 and 
5 millirems per hour at a distance of 10 feet (3 m) from the source. It was assumed that 
the remaining component weights would be unimportant with respect to these magnitudes. 

The temperature range in figure 11 was extended beyond the range of initial interest 
in order to determine potential benefits. Both fuel -block and radiator weight decrease 
as turbine inlet temperature is increased in the range from 1800° to 2200° R (1000° to 
1222° K). However, the fuel block with its assumed 2460° R (1367° K) temperature 
limitation requires increasingly large surface areas (see appendix D). This results in 
increasing shield weights about 1900° R (1056° K). Hence, the total of these component 
weights exhibits minimums for dose rates of 5 and 1 millirem per hour, respectively, 
at approximately 2050° and 2000° R (1140° and 1110° K). 

Compressor inlet temperature. - Figure 12 presents the effect of changing radiator 
area and compressor inlet temperature on component weights. The lower set of curves 
shows that as compressor inlet temperature (radiator outlet temperature) is increased, 
radiator weight for a given area rapidly decreases and levels off at a given value. Shield 
weight (for a dose rate of 5 millirems per hour at 10 feet (3 m)), fuel -block weights, and 
the sum of the three component weights are presented in the top part of the figure. 
Decreasing compressor inlet temperature increases cycle efficiency at the expense of 
increasing radiator area (see Thermodynamic Variations). At each compressor inlet 
temperature, there is a minimum value of radiator area (i. e. , the prime radiator area) 
required to reject the cycle heat. From a weight standpoint, a slightly larger radiator 

would be much lighter. For example, at 540° R (300° K), the minimum-area^ radiator is 

2 

approximately 430 square feet (40 m ) and 1600 pounds (720 kg). An increase in radiator 

2 

The drop between the bulk liquid temperature and the radiating surface temperature 
is included in the radiator weight calculations. Hence, the "exact" minimum or prime 
area is approached as radiator weight becomes large. 
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2 

area to 450 square feet (41. 8 m ) decreases radiator weight to 800 pounds (360 kg). For 
some missions, existing vehicle structure might be used as a part of the radiator struc- 
ture, and, hence, area itself might be more restrictive than weight. As compressor 
inlet temperature is decreased, the minimum required radiator area rapidly increases. 
Selection of a low value for compressor inlet temperature could preclude use of the power 
system for some surface -area-limited applications. 

The minimum total weight shown in figure 12 for 1000 square feet (92. 9 m^) of 
radiator area is 3400 pounds (1540 kg) and occurs at a compressor inlet temperature of 
approximately 495° R (273° K) or 45° above the assumed sink temperature. As surface 


. 09r 


,06b 


• 03r 




Figure 13. - Variations of heat-exchanger weights and volumes with 
recuperator effectiveness. Gross shaft power, 6.875 kilowatts. 

{All other initial reference conditions of table I apply. ) 
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area is reduced, the compressor inlet temperature for minimum total weight increases; 
for example, in decreasing from 1000 to 800 to 600 to 400 square feet (92. 9 to 74. 3 to 
55. 8 to 37.2 m^), the compressor inlet temperatures were 495°, 500°, 515°, and 565° R 
(273°, 278°, 286°, and 314° K), respectively. 

Recuper ator effectiveness . - Figure 13 shows the variations in weights and volumes 
of recuperator, heat source, and heat-sink heat exchanger over a range of recuperator 
effectiveness (0. 87 to 0.95). In addition, radiator weight for 650 square feet (60. 3 m^) 
of radiator area is shown. Pressure-drop distributions were determined in appendix E 
for minimum combined heat-source, recuperator, and heat-sink heat -exchanger weights. 

Increasing recuperator effectiveness increases cycle efficiency and decreases radi- 
ator weight for a fixed area but increases heat- exchanger weight and volume. For the 
range in effectiveness from 0. 875 to 0. 950, recuperator weight increases from approxi- 
mately 82 to 170 pounds (37 to 77 kg), while recuperator core volume increases from 



Recuperator effectiveness, E r 


Figure 14. - Effect of recuperator effectiveness on major component weights. 
Gross shaft power, 6.875 kilowatts. (All other initial reference conditions 
of table I apply. ) 
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o 

approximately 0. 52 to 1. 52 cubic feet (0. 0147 to 0. 043 m ). Increasing effectiveness 

has relatively little effect on the heat-source and heat-sink heat -exchanger weights and 

volumes because of the opposing effects of reduced heat loads and decreased pressure 

drops. The total weight for these components shows a relatively small variation over 

this range of recuperator effectiveness, but with an insensitive minimum of 489 pounds 

(222 kg) at an effectiveness of approximately 0. 90. The total combined heat-source, 

recuperator, and heat-sink heat-exchanger volumes continually increase from about 

3 3 

1. 5 cubic feet (0. 0424 m ) at an effectiveness of 0. 875 to 2. 68 cubic feet (0. 0759 m ) at 

0. 950. In spite of the relatively large percentage recuperator weight changes, these 
variations are small compared with the fuel -block weights. 

Figure 14 presents variations of the major component weights and their sums with 
recuperator effectiveness. Because of the preponderance of heat-source weights, the 
range of effectiveness was extended up to 0. 990 in order to establish minimum weight. 
Weights of fuel block, shield (for dose rates of 1 and 5 mrem/hr), combined heat 
exchanger, and the corresponding total of major component weights are presented along 
with the same total minus shield weights. The total unshielded weight has a minimum 
at a recuperator effectiveness of about 0. 95. However, the unshielded weight does not 
vary appreciably for values of effectiveness from 0. 90 to 0. 98. The total weights for 
dose rates of 1 and 5 millirems per hour showed minimums at values of recuperator ef- 
effectiveness of approximately 0.97 and 0.96, respectively. 

System loss pressu re ratio . - Effects of variation in system loss pressure ratio of 
0. 90 to 0. 95 on heat-exchanger weights and volumes are presented in figure 15. Recuper- 
ator, heat-source and heat-sink heat-exchanger weights and volumes, as well as radiator 

2 

weight for 650 square feet (60. 3 m ) of area, and the total weights and volumes are all 
shown. The pressure-drop distributions as a function of total system pressure loss are 
presented in appendix E. Increasing loss pressure ratio increases cycle efficiency and hence 
reduces fuel-block and radiator weights but at the expense of heat -exchanger weights and 
volumes. Furthermore, small values of total system pressure loss require small compo- 
nent pressure drops and, in general, require more stringent designs. 

Each of the three heat exchangers shows relatively small increases in weight with 

2 

increasing loss pressure ratio. Addition of the radiator weight (for an area of 650 ft 

o 

or 60. 3 m ) to the sum of the other heat-exchanger weights results in a minimum com- 
bined weight of about 430 pounds (195 kg) at a system loss pressure ratio of 0. 925. 

Again, these variations are small compared with the heat -source weights. 

Both the recuperator and heat-sink heat exchanger show relatively small increases 
in volume (20 percent) with increasing loss pressure ratio, while the heat-source heat 
exchanger shows nearly a threefold increase in volume from a loss pressure ratio of 
0. 900 to 0. 950. For minimum combined weight (see appendix E), the percentage pres- 
sure drop allocated to the heat-source heat exchanger varies from 3. 7 percent at a loss 



Loss pressure ratio, L 


Figure 15. - Effect of loss pressure ratio on 
heat-exchanger weights and volumes. 
Gross shaft power, 6. 875 kilowatts. (All 
other initial reference conditions of table I 
apply. ) 
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Figure 16. - Effect of loss pressure ratio on 
major component weights. Gross shaft 
power, 6.875 kilowatts. (All other initial 
reference conditions of table I apply. ) 


pressure ratio of 0. 90 to 1. 6 percent at a loss pressure ratio of 0. 95. The total heat- 

3 

exchanger volume continuously increases from about 1. 7 cubic feet (0. 048 m ) at a loss 

3 

pressure ratio of 0. 90 to about 3 cubic feet (0. 1 m ) at 0. 95. 

Effects of variation in loss pressure ratio from 0. 90 to 0. 95 on the major compo- 
nent weights and their sums are presented in figure 16. The total heat-exchanger weight 
from figure 15, fuel-block weight, unshielded weight, a band of shield weights, and a 
band of the total weights are shown in figure 16. The sums show continously decreasing 
weight with increasing loss pressure ratio. 


Selection and Performance Estimates 


Results of the parametric variations are now summarized with a selection of Brayton- 
cycle parameters suitable for a radioisotope power module to be designed in the electric 
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power range from 2 to 10 kilowatts. This is followed by performance estimates for the 
conceptual powerplant at 2, 6, and 10 kilowatts. 

Selections . - The Brayton-cycle parameter selections and major reasons for the 
selections were as follows: 

(1) Working -fluid molecular weight, 83.8: for the same pressure levels, this molecular 

weight compared with 39. 9 resulted in a larger diameter turbine and compressor 
at lower rotational speeds and potentially higher cycle efficiency (fig. 10) 

(2) Rotational speed, 36 000 rpm (3770 rad/sec): a reasonable compromise among 

pressure level, turbine and compressor diameter, and cycle efficiency (fig. 6) 

(3) Turbine inlet temperature, 2060° R (1144° K): this value resulted in high cycle 

efficiency, low minimum radiator area, and near -minimum total system weight 
(figs. 4 and 11) 

(4) Compressor inlet temperature, 540° R (300° K): a reasonable compromise 

which resulted in relatively low minimum radiator area without large penalties 
in total system weight (fig. 12) 

(5) Recuperator effectiveness, 0. 95: resulted in near -minimum total system weight 

and still yielded relatively compact heat-exchanger designs (figs. 13 and 14) 


p a, net ' 6.0 kW 
^sy-0.23 
f - 1200 Hz 


1711° R 
14.41 psia 


Controls: 
0.40 kW 


2060° R 
25.28 psia 


Heat source: 

Qi =26.2kW 



Helium-xenon: 
• 83.8 


m = 0.963 Ib/sec 


1663° R 
25.93 psia 
1663° R 
26.01 psia 


540° R 
13.93 psia 


1711° R 
14. 32 psia 


T — — * 

[jSfEj- 0,95^1 _ 

747° R 

! 111., ail 527° R ' 



CD-9635 


(a) U. S. Customary Units. 

Figure 17. - Schematic diagram of conceptual powerplant at 6 kilowatts net power. 
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(6) Loss pressure ratio, 0. 92: relatively low system weight and relative ease of 

heat-exchanger component designs (fig. 16) 

(7) Compressor pressure ratio, 1. 90: for the selected recuperator effectiveness 

and loss pressure ratio, this value should result in near -maximum efficiency 
and reasonable pressure levels (see figs. 5, 7, and 8 for trends) 

(8) Compressor specific speed, 0.7750: with the preceding selections, this value 

should result in high efficiency and comparatively high bearing cavity pressure 
for the 2 -kilowatt power level but not prohibitively high with respect to generator 
windage loss at the 10-kilowatt power level (appendix F) 

Performance. - The preceding selections result in a cycle thermal efficiency of 
0. 360 and a specific compressor inlet pressure of 1. 62 psi per kilowatt (11. 2 (N/m )/kW). 
Turbine- and compressor -tip diameters were 4. 9 and 3. 8 inches (12. 42 and 9. 65 cm), 
respectively, while the corresponding efficiencies were predicted to be 0. 85 and 0. 76. 

A schematic diagram of the conceptual powerplant employing the parametric selec- 
tions is shown in figure 17 for a net unconditioned power output of 6 kilowatts. The 
alternator was assumed to be a four-pole Lundell type with a smooth 3. 6 -inch - 
(9. 14-cm-) diameter rotor (see appendix F for a comparison to homopolar -inductor -type 


= 6.0 kW 


950° K 


95.9 kN/m 2 

/-Helium-xenon: 

’ A 

m = 0.436 kg/sec 



CD-9635 


(b) SI Units. 
Figure 17. - Concluded. 
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alternators). Resulting output frequency is 1200 hertz. The alternator was sized for 
10 kilowatts of net system output and was estimated to have an electromagnetic effi- 
ciency of 0. 89 at the 6 -kilowatt power level. A gas mass flow rate of 0. 963 pound per 
second (0. 436 kg/sec) produces a gross shaft power of 8. 56 kilowatts. Estimates of 
alternator windage and bearing friction and an allowance for shaft-seal losses amounted 
to 1. 38 kilowatts; net useful shaft power is then 7. 18 kilowatts. Gross alternator output 
is 6. 4 kilowatts. An allowance of 400 watts was made for the combination of speed con- 
trol, voltage regulator and exciter, parasitic load control, and liquid-coolant pump-motor 
power. 

Required cycle thermal input is 23. 8 kilowatts. An allowance of 1. 2 kilowatts was 
made for system thermal losses, resulting in a fuel block thermal power requirement of 
25. 0 kilowatts at the end of the assumed 5-year mission. Allowing for isotope decay 
over the 5 years and 90 days (assumed launch-delay time) results in an initial thermal 
inventory requirement of 26.2 kilowatts. Hence, initial overall system efficiency at the 
6 -kilowatt net power level is 0. 23. 

The parametric selections result in a compressor inlet pressure of 13. 89 psia 
(95. 7 kN/m ) for the 6-kilowatt net unconditioned power level. A compressor pressure 
ratio of 1. 90 and loss pressure ratio of 0.92 results in the turbine pressure ratio of 
1. 748. Pressure drop distributions among the heat -exchanger components are as follows 
(appendix E): recuperator, 3. 5 percent; heat-source heat exchanger, 2. 5 percent; 
piping, 1. 5 percent; and heat-sink heat exchanger, 0. 5 percent. With the assumed heat- 
sink heat-exchanger effectiveness of 0. 95 and capacity-rate ratio of 0. 87, the minimum 

specific radiator area required was approximately 65 square feet per kilowatt electric or 

2 2 

6. 04 square meters per kilowatt electric (390 ft , or 36.2 m , for 6-kW net power). 

Table III presents breakdowns of the system efficiency calculations at net system 
power levels of 2, 6, and 10 kilowatts. The same combined rotating unit was to be used 
at each power level, but new heat -transfer components were employed to maintain the 
design parameters selected in the preceding section. The alternator was sized to pro- 
duce up to 10 kilowatts of net power, while the journal bearings were sized at the 
2 -kilowatt power level (lowest cavity pressure). Changes in compressor and turbine 
efficiencies over this power range from Reynolds number effects were estimated to be 
small and were neglected. Hence, the compressor and turbine pressure ratios and effi- 
ciencies were assumed to be constant. 

Among the parasitic losses, alternator windage dominated at the 6- and 10-kilowatt 
power levels, while the bearing friction was the largest of these losses at the 2 -kilowatt 
level. For constant size and Reynolds number, alternator windage varies directly with 
gas pressure. The sum of the journal- and thrust -bearing frictions was estimated to be 
fairly constant over the power range. Seal leakage was accommodated by assuming that 
the compressor mass flow rate was 1 percent greater than the turbine mass flow rate. 


33 



TABLE m. - SYSTEM EFFICIENCY ESTIMATES 



Power output, 



kW 



2 

6 

10 

Cycle parameters: 




Initial fuel block power, kW 

10.9 

26.2 

41. 5 

Thermal loss, kW 

.5 

1.2 

1.9 

Cycle input power, kW 

9.9 

23. 8 

37.7 

Cycle thermal efficiency 

.360 

. 360 

. 360 

Gross shaft power, kW 

3. 56 

8. 56 

13. 57 

Parasitic losses: 




Alternator windage, kW 

.35 

.64 

. 88 

Bearing friction, kW 

. 52 

. 53 

. 54 

Seal leakage, kW 

. 09 

.21 

. 34 

Net shaft power, kW 

1 2.60 

7. 18 

11. 81 

System parameters: 




Alternator efficiency 

.83 

.89 

.90 

Gross alternator power, kW 

2. 15 

6. 40 

10.64 

Radiator coolant pump power, kW 

. 06 

. 14 

.22 

Control power, kW 

. 09 

.26 

. 42 

Net system power (unconditioned), 

2. 00 

6. 00 

10. 00 

kW 




System efficiency 

0. 183 

0. 229 

0. 241 


Alternator design point (10 kW) electromagnetic efficiency was estimated at 0. 90, 
with a degradation to 0. 89 at 6 kilowatts and 0. 83 at 2 kilowatts. The remaining allow- 
ances were for the power requirements of the system itself. System efficiency was 
estimated to be 0. 18, 0. 23, and 0. 24 at net system power levels of 2, 6, and 10 kilowatts, 
respectively. 

It appears that tailoring the combined rotating unit for a specific power level rather 
than the range from 2 to 10 kilowatts would result in relatively small efficiency improve- 
ments. If the design were only for the 10-kilowatt level, it was estimated that about 
180 watts could be saved by reducing bearing size. Hence, the sizing for the 2-kilowatt 
level resulted in a l| percent penalty in system efficiency at the 10-kilowatt level. 
Conversely, if the design were only for the 2-kilowatt level, it was estimated that about 
0. 8 improvement in alternator electromagnetic efficiency could be obtained. Hence, the 
sizing for the 10-kilowatt level resulted in a 7. 8 percent penalty in system efficiency at 
the 2 -kilowatt level. 

Table IV breaks down the system weight at the net system power levels of 2, 6, and 
10 kilowatts. On the basis of minimum combined fuel -block and shield weight, the fuel 
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TABLE IV. - SYSTEM WEIGHT ESTIMATES 


(a) U. S. Customary Units 


Power output, 
kW 



2 

6 

10 

Heat source: 




Fuel block, lb 

527 

1226 

2008 

Reentry body, lb 

150 

345 

565 

Emergency heat sink, lb 

65 

156 

248 

Subtotals: 




lb 

742 

1727 

2821 

lb/kW 

371 

288 

282 

Power conversion unit: 




Combined rotating unit, lb 

125 

125 

125 

Source heat exchanger, lb , 

122 

91 

107 

Recuperator, lb I 

260 

210 

240 

Sink heat exchanger, lb 

123 

55 

60 

Ducting, lb 

90 

80 

85 

Structure and insulation, lb 

200 

160 

180 

Controls, lb 

175 

220 

270 

Subtotals: 




lb 

1095 

941 

1067 

lb/kW 

548 

157 

107 

Radiator: 




Area, ft^ 

180 

540 

900 

Weight: 




lb 

150 

282 

! 495 

lb/kW 

75 

47 

50 

Unshielded totals: 




lb 

1987 

2950 

4383 

lb/kW 

994 

492 

438 

Shield (5 mrem/hr): 




lb 

1190 

2150 

3050 

lb/kW 

595 

360 

305 

Totals: 




lb 

3177 

5100 

7433 

lb/kW 

1588 

850 

743 



TABLE IV. - Concluded. SYSTEM MASS ESTIMATES 


(b) SI Units 



Power output, 
kW 


2 

6 

10 

Heat source: 




Fuel block, kg 

240 

556 

912 

Reentry body, kg 

68 

157 

256 

Emergency heat sink, kg 

30 

71 

113 




Subtotals: 




kg 

338 

784 

1281 

kg/kW 

169 

131 

128 

Power conversion unit: 




Combined rotating unit, kg 

56.7 

56.7 

56.7 

Source heat exchanger, kg 

55.4 

41. 3 

48.6 

Recuperator, kg 

118 

95. 3 

109 

Sink heat exchanger, kg 

55.8 

25. 0 

27.2 

Ducting, kg 

40.9 

36. 3 

38. 6 

Structure and insulation, kg 

90.8 

72. 6 

81.6 

Controls, kg 

79. 5 

100 

123 

Subtotals: 




kg 

497 

427 

485 

kg/kW 

249 

71.2 

48.5 

Radiator : 




Area, m^ 

16.9 

50.6 

84. 5 

Weight: 




kg 

68. 1 

128 

225 

kg/kW 

34.0 

21. 4 

22. 5 

Unshielded totals: 




kg 

903 

1340 

1990 

kg/kW 

451 

224 

199 

Shield (5 mrem/hr): 




kg 

540 

976 

1385 

kg/kW 

270 

163 

139 

Totals: 




kg 

1440 

2320 

3370 

kg/kW 

720 

386 

337 


block was assumed to have one row of capsules with each capsule void volume sized for 
a maximum internal pressure of 12 500 psi, or 86. 1 meganewtons per square meter 
(appendix D). As shown in table IV, the total heat-source weight increased with 
increasing system power level, while the specific weight decreased. These specific 
weights represent from about 1/3 to 2/3 of the total unshielded system specific weights 
between 2 and 10 kilowatts, respectively. 

The combined rotating unit was estimated to weigh 125 pounds (56. 8 kg). Each of the 
three heat exchangers was sized for its respective heat loads, pressure levels, and 
percent pressure drops (all three of which varied with power level; see appendix E). 

Since the heat-source heat -exchanger weight was small with respect to total system 
weight, this component was sized for minimum volume rather than minimum weight. 
Presumably, the relatively small penalty (approximately 30 lb, or 13. 6 kg) in the 
heat -source heat -exchanger weight would be more than compensated for by shield weight 
reduction and ease of spacing arrangements within the shielded volume. 

The total weights for the power conversion unit were nearly constant for the three 
power levels. Opposing effects of the variations in heat load and pressure levels resulted 
in no major variation in heat -exchanger and hence total power conversion unit weight with 
power level. Specific weights, on the other hand, decreased by more than 70 percent in 
going from 2 to 6 kilowatts, and by an additional 32 percent in going from 6 to 
10 kilowatts. The power conversion unit specific weight ranged from about 1/2 to 1/4 
of the total unshielded system specific weight between 2 and 10 kilowatts, respectively. 

A radiator area of 90 square feet per kilowatt electric (8. 36 m /kWe) was assumed 

so that the radiator weights in table IV are near minimum (see appendix E). Minimum 

specific area requirements were approximately 65 square feet per kilowatt electric 

2 

(6. 04 m /kWe). Structural aspects of the radiators were not considered. 

The shield weight estimates were based on a dose rate of 5 millirems per hour at a 
distance of 10 feet (3 m) from the isotope heat source. The absolute weight of this com- 
ponent increased with increasing system power output. The shield specific weights 
represented approximately 40 percent of the total shielded system specific weight. Both 
shielded and unshielded system specific weights decrease by more than a factor of 2 over 
the power range from 2 to 10 kilowatts. 

A conceptual arrangement of the power conversion unit is shown in figure 18 for a 
net power output of 6 kilowatts electric. The overall package size is about 50 by 40 by 
20 inches (125 by 100 by 50 cm). The hot gas from the heat-source heat exchanger 
enters the turbine scroll, leaves axially, and is then turned into the triangular end sec- 
tion of the recuperator. The flow progresses through the recuperator and heat -sink heat 
exchanger and is then turned back to an axial direction before entering the compressor. 
The compressed gases leave the compressor scroll and are ducted through the alternate 
passages of the recuperator and are fed back to the heat-source heat exchanger. 
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Figure 18. - Conceptual 6-kilowatt power conversion unit arrangement. 

The core dimensions of the recuperator were 12. 3 by 12. 3 by 19. 6 inches (31. 2 by 
31. 2 by 49. 8 cm) while the heat-sink heat exchanger had the dimensions 12. 3 by 8.7 by 
9. 0 inches (31. 2 by 22. 1 by 22. 9 cm). The combined rotating package was estimated to 
have a cylindrical envelope, 23 inches (58. 4 cm) long and 13 inches (33 cm) in diameter. 

The fuel block (see fig. 21, appendix D) for the 6-kilowatt-electric net power level 
required a face area of 7 square feet, or 0. 65 square meter (31. 8 by 31. 8 in. , or 80. 8 
by 80. 8 cm) and was 3. 6 inches (8. 9 cm) thick. One row of capsules, each with an inner 
diameter of 2. 40 inches (6. 10 cm) was used. The heat-source heat exchanger had the 
dimensions 31. 8 by 31. 8 by 0. 9 inch (80. 8 by 80. 8 by 2. 29 cm). 

SUMMARY OF RESULTS 

Parametric analysis of a conceptual plutonium 238, Brayton-cycle space powerplant 
yielded the following results: 
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1. The following major Brayton-cycle design variables were selected for a system 
to produce rated power from 2 to 10 kilowatts: 


Turbine inlet temperature, °R (°K) 2060 (1144) 

Compressor inlet temperature, °R (°K) 540 (300) 

Working fluid molecular weight 83. 8 

Shaft rotational speed, rpm (rad/sec) 36 000 (3770) 

Compressor pressure ratio 1. 90 

Compressor specific speed 0. 775 

System loss pressure ratio 0. 92 

Recuperator effectiveness : 0. 95 


2. One combined rotating unit was sized to operate over the range of rated powers 
from 2 to 10 kilowatts, while the gas pressure level was varied. Performance penalties 
for use of one combined rotating unit over this power range (compared with redesign for 
each power) were estimated to be small. This rotating unit consists of single-stage 
radial -flow turbine and compressor and a smooth-rotor four -pole alternator, all sup- 
ported by hydrodynamic gas bearings. Resulting tip diameters for the turbine, alterna- 
tor, and compressor were 4. 9, 3. 6, and 3. 8 inches (12. 4, 9. 1, and 9. 7 cm), respec- 
tively. Efficiencies for the turbine and compressor were predicted to be 0. 85 and 0. 76, 
respectively, while variations due to power level (or Reynolds number) were estimated 
to be small and were neglected. 

3. Performance was estimated for 2, 6, and 10 kilowatts of net power output. The 
heat -exchanger components were resized at each power level. For this range of power, 
cycle thermal efficiency was 0. 36, and the ratio of compressor inlet pressure to gross 
shaft power was 1.62 psi-kilowatt (11. 18 (kN/m 2 )-kW). 

4. For minimum combined weight of heat source and nuclear shielding, radioisotope 
capsule void volume was sized for a maximum internal helium pressure of 12 500 psi 
(86. 1 MN/m 2 ). 

5. Overall system efficiency (the ratio of net electrical output to initial isotope 
thermal power), unshielded system specific weights, and shielded specific weights were 
estimated as follows: 


Net 

power, 

kW 

Overall 

efficiency 

Unshielded 
specific weight 

Shielded (5 mrem/hr) 
specific weight 

lb/kWe 

kg/kWe 

lb/kWe 

kg/kWe 

2 

0. 18 

990 

449 

1590 

724 

6 

.23 

490 

222 

850 

386 

10 

.24 

440 

200 

740 

336 
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o 

These weights assumed 90 square feet per kilowatt electric (8. 36 m /kWe) of radiator 
area and resulted in near -minimum radiator weight. 

6. Minimum radiator area was approximately 65 square feet per kilowatt electric 
(6. 04 m 2 /kWe). 

7. An envelope of the power conversion unit, excluding the heat-source heat exchan- 
ger, was estimated at the 6-kilowatt power level to be 20 by 40 by 50 inches (50 by 100 
by 125 cm). 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, January 19, 1968, 

120-27-03-28-22. 
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APPENDIX A 


SYMBOLS 



2 2 

prime radiator area, ft ; m 
frontal area of fuel block, ft^; 
constant in eq. (B7) 


constant in eq. (B7) 

specific heat at constant pressure, Btu/(lb mass)(°R); J/ (kg)(°K) 


capacity rate ratio of heat-sink heat exchanger (ref. 


12), mCp/(mCp)^ 


axial velocity component, ft/ sec; m/sec 


3 D [gJ(Ah) id ] 1/4 

specific diameter , — = — , dimensionless; 

12 VQ 



D t tip diameter, in. ; cm 

E r recuperator effectiveness 

E g heat-sink heat-exchanger effectiveness 

f frequency, Hz 

p p 

g acceleration due to gravity, 32. 2 ft/ (sec ;(lb mass/lb force); 9. 81 m/ (sec )(kg/N) 

h gas enthalpy, Btu/lbmass; J/kg 

J mechanical equivalent of heat, 778 (ft)(lbforce)/Btu; (1 Nm/J) 

constant in eq. (B9), (in.) - *; (cm) - * 

K m constant in eq. (B8) 

L system loss pressure ratio, (p^/pg .st>(P4 /p 5 ) 

o 

L turbine exit loss parameter, (c 2 /V.) 

e ^ J 

L. heat-source insulation loss, kW 

Lj journal bearing friction loss, kW 

L shaft seal leakage loss, kW 

s 

Lj. thrust bearing friction loss, kW 

3 For common usage (refs. 13 and 14, e.g. ) multiply dimensionless D g by 0. 420 to 

obtain (sec*/'Vft*/ 3 4 )(lbforce/lbmass)*/ 4 . 
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w 

M* 


M 


w 


m 


N 

N 


s 


n 

P 

P sh 

p 

sn, net 
P 

Ap/p 

Q 

Q. 

Tin 

R 

Re 

r 

T 

T s 

U t 

v ; 

Q! 

/3 

y 

A 

e 


alternator windage loss, kW 
characteristic tip Mach number 
gas molecular weight 
gas mass flow rate, lb mass/hr; kg/hr 
shaft rotational speed, rpm; rad/sec 


specific speed 


4 


ttnVq 

30[gJ(Ah). d ] 3/4 


, dimensionless; 



exponent in eq. (B7) 
power, kW 

gross shaft power, kW 
net shaft power, kW 

O 

pressure, psi; N/ cm 
component pressure drop, percent 
volume flow rate, ft / sec; m / sec 
cycle thermal input, kW 

gas constant, (ft)(lbforce)/(°R)(lbmass); J/(kg)(°K) 
Reynolds number 
machine pressure ratio 
temperature, °R; °K 

effective radiator sink temperature, 450° R; 250° K 
tip speed, ft/sec; m/sec 
ideal jet speed, ft/sec; m/sec 

pressure drop distribution ratio, (Ap/ P^pg/ (Ap/p) lpg 

p 

compressor slip factor, gJ(Ah) /Ur 

L Cr 

ratio of specific heats 

difference operator 

surface hemispherical emittance 


For common usage (refs. 13 and 14, e.g.) multiply dimensionless N by 129 to 
obtain |(rpm)(ft) 3 / 4 /(sec) 1//2 J (lbmass/lbforce) 3//4 . 
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77 efficiency 

6 area multiplication factor for radiation interchange 

p gas viscosity, lb mass/(ft)(sec); kg/(m)(sec) 

v turbine blade- to ideal -jet speed ratio 

O O 

p gas density, lb mass/ft ; kg/m 

a Stefan-Boltzmann constant, 0. 173x10”® Btu/(hr)(ft 2 )(°R) 4 ; 5. 67x10”® 

j/(sec)(m 2 )(°K) 4 

Subscripts: 
a alternator 

c compressor 

cr critical 

cy cycle 

em electromagnetic 

g gross 

hps high pressure side 

id ideal 

k controls 

l liquid 

lps low pressure side 

max maximum 

net net 

0 initial 

opt optimum 

p pump 

st static 

sy system 

t turbine 

1 turbine inlet condition 

2 turbine exit condition 

3 heat-sink heat -exchanger (gas side) inlet condition 


43 



4 compressor inlet condition 

5 compressor exit condition 

6 heat-source heat-exchanger inlet condition 
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APPENDIX B 


RADIAL TURBOMACHINERY EFFICIENCY MODELS 

By dimensional analysis, reference 15 shows that only four parameters are required 
to describe completely the performance of compressible flow turbomachinery, provided 
that complete physical similarity can be maintained. Quoting from reference 15, 
"Complete physical similitude implies (1) geometrical similitude, which means that the 
linear dimension ratios are everywhere the same, that is, shapes are similar regardless 
of size; (2) kinematic similitude, which means that velocity ratios are the same, i. e. , 
the velocity triangles representing flow conditions are similar; and (3) dynamic 
similitude, which means that the ratios of the different forces are everywhere the same. 
For fluid flow, the streamlines are everywhere similar. " In practice, complete simili- 
tude is difficult to attain. This is especially true for small machinery where manufac- 
turing techniques necessary to attain similar clearance ratios are not practical. (These 
techniques could be achieved but currently do not warrant the cost. ) 

The empirical models of this analysis were based on the four similarity parameters 
of references 13 and 14 with minor modifications and the inclusion of a "size effect. " 
Hence, the total efficiencies of both turbine and compressor were assumed to be functions 
only of specific speed, specific diameter, Reynolds number, characteristic tip Mach 
number, and tip diameter, or expressed functionally, 

V = f (N g , D g , Re, M*, D t ) (Bl) 


The model assumed that the compilation of best efficiency points for radial -flow 
turbines and compressors as a function of specific speed, presented in reference 16, 
represented the maximum obtainable efficiency, that is, the efficiency obtainable when 
the remaining parameters were at their optimum values or beyond their range of influ- 
ence. Hence, from reference 16, 


r ? = %iax =f < N s ) 


(B2) 


for 


Re > (Re) cr 

M* < (M*) 

x 'cr 
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D t - cv 


cr 


D 


s 



) 

opt 


It was assumed that the optimum specific diameter would result from the specification of 
rotational speed and a blade- to ideal -jet speed ratio of 0. 70 for the turbine, and rota- 
tional speed and a 0. 85 value of slip factor for the straight radial -bladed compressor. 
Symbolically, this is represented as 


where for a turbine, 


<D S ) =f(D t ) 
opt 


(B3) 


720 U, / 200 U.N 


7rN 


N 


(B4) 


U, 




= 0.70 y2gJ(Ah) t id ^ 


st 


(B5) 


and for a compressor, 


U, 


t, c 




(Ah), 


i 


gJ_ 

0. 85 


(Ah) 


(B6) 


A critical Reynolds number of 10 was assumed and efficiency corrections were 
assumed to follow the general form (ref. 15) 

/ F7\ ri"1 


^ = 1 - 


a + b 


f 10 

[Re. 


(B7) 


for Re < 10. Mach number (or pressure ratio, since the two are related) effects were 
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assumed to be linear above critical values and have the general form 


A,, = - K m (M* - M* r ) (B8) 

for M* > M* A critical tip diameter of 6 inches (15 cm) was assumed with efficiency 

V/l 

corrections empirically fitted to 


At, = K d (6. 0 - D t ); [at, = - K d (15. 24 - D t )] (B9) 

for D t < 6. 0 p t < 15.24). 

From the results of the 6-inch (15-cm) turbine and compressor investigations 
(ref. 6), the following coefficients were determined for equation (B7) (the effect of 
Reynolds number). For the turbine, 

a = 0. 3 b = 0.7 n = 0. 20 

and for the compressor, 

a = 0, b = 1.0, n .= 0. 06 

Comparisons among analytical and experimental results of references 5, 13, 14, 
and 17 led to selection of the constants of the Mach number correction (eq. (B8)) for the 
compressor. These were K m = 0. 075 and M£ r = 0. 5. For the Mach number effect on 
turbines, the analytical results of reference 17 were employed. The constants were 
K m = 0. 070 and M* r = 0. 25 for the range of 0. 25 < M* < 0. 75, and K m = 0. 320 and 
M cr = ^ f° r ^ ^5- Using the preceding corrections and the peak performance 

results shown in table V, the remaining efficiency differences were allocated to the 
"size effect. " The resulting coefficients were 

K d = 0. 0125; [K d = 0. 00403] 

for compressors with D^_ < 6. 0 (D^. < 15. 24) and 

K d = 0. 0085; [x d = 0. 00335] 

for turbines with > 6. 0 (D^. > 15. 24). 

The procedure then in calculating a turbine or compressor efficiency was to assume 
a rotational speed and a specific speed and to determine the maximum efficiency from the 
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empirical prediction of reference 16 (eq. (B2)). Then, this calculated efficiency was 
corrected for Reynolds number (eq. (B7)), Mach number (eq. (B8)), and size (eqs. (B4) 
and (B9)) in that order. 


TABLE V. - COMPARISON OF RADIAL-FLOW 
TURBOMACHINERY AT PEAK EFFICIENCY 


Parameter 

Turbines 

Compressors 

Tip diameter, D,; 





. 2 1 
in. 

6. 02 

4. 59 

6. 0 

3.2 

cm 

15.3 

11.7 

15.2 

8. 12 

Total pressure ratio 

1. 565 

1. 555 

2.33 

2. 00 

Static pressure ratio 

1.610 

1. 605 



Total efficiency 

0. 88 

0.864 

0.745 

a 0. 750 

Static efficiency 

0. 83 

0. 810 



Specific speed, N g 

0. 688 

0.734 

0.631 

0. 691 

Specific diameter, D g 

2.95 

2.78 

3. 89 

3.46 

Reynolds number, Re 

' 3. 5xl0 5 

4. 4xl0 5 

1. 39xl0 6 

3. 12xl0 6 

Tip Mach number, M* 

0. 582 

0.579 

1. 104 

0.956 

Exit loss parameter, L g 

0. 059 

0. 064 



Slip factor, jB 



0. 89 

0.93 


a Estimated as 0. 05 lower than measured thermal efficiency. 
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APPENDIX C 


CYCLE CALCULATION PROCEDURE 

A digital machine computer code was constructed to analyze the thermodynamics of 
the Brayton cycle system utilizing the efficiency models presented in appendix B. The 
equations and procedure of this analysis are presented. 

The input parameters basically were those listed in table I. For gas properties, 
the analysis assumed perfect gas relations. Compressor specific work was calculated 
from 


(ah) = fell 

^c 


'£5' 

V P 4 


(y-l)/y 


- 1 


(Cl) 


where 


y - 1 J \ 3 


(C2) 


and 


R 


1545 


M 


w 


R = 


8317 




w. 


(C3) 


Initially, an arbitrary value was assumed for rj , while and p^/p^ were specified 
in the input. The pressure ratio available to the turbine was given by 


P2. st \P4/ 


(C4) 


and the turbine specific work was determined from 


< Ah >t = \ st c p T l 


1 - 


Pi 


-(y-l)/r 


V P 2, st/ 


(C5) 
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Initially, efficiency was an assumed value. For a fixed gross shaft power, the cycle 
mass flow rate was given by 


M = 


3415P sh 
(Ah) t - (Ah) c ’ 


M = 


1550P sh 
(Ah) t - (Ah) c 


(C6) 


The next steps sized the turbine and compressor and computed the component effi- 
ciency by iterating until the assumed efficiencies agreed with the calculated efficiencies 
within ±0. 0005. Compressor tip speed was calculated from 


and tip diameter from 




720U, „ 
t,e 

7rN 



100N 
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With rotational speed and compressor specific speed assigned, the compressor char- 
acteristic volume flow (inlet) was obtained from 



Hence, the inlet density and pressure were evaluated as 


(C9) 


P 4 = 


m 

3600Q 4 


(CIO) 


and 


P 4 = 


4 

144 


[p 4 = P4 RT 4] 


(Cll) 


The characteristic tip Mach number (ratio of tip speed to the critical speed of sound at 
the inlet) and Reynolds number were obtained from 
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(C12) 


Re = P A £V : 
C 12M 4 


PdPt c D t c 
Re = _ h c I± S . 

c 100u 4 


(CIS) 


Fluid viscosities (as a function of temperature) and maximum component efficiencies (as 
a function of specific speed) were handled as tables which were interpolated in appro- 
priate subroutines. With the preceding calculations used as input, an efficiency for the 
compressor was calculated from equation (B7). 

Turbine tip speed and diameter were calculated from 




(C14) 


72 OIL 


t ’ t 100N 


(C 15) 


where 


v V 2qJ(Ah) t,id,st 


(C 16) 


; t, id, st 


\st 


(C17) 


The turbine characteristic volume flow Qg ^ (exit static conditions) was developed from 
the following relations: 


T 2 = T 1 - 


(CIS) 
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(C19) 
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where a = (Ap/p) hps /(Ap/p) lpg 
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The parameters L, L g , and a were input quantities. 

Turbine total efficiency was related to its static efficiency by 

\st 
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Turbine specific speed was calculated from 
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The characteristic blade tip Mach number and Reynolds number based on inlet conditions 
were evaluated from 


M t* = 



(C26) 


and 


Re, 


p l U t,t D t,t, 

12 Mi 


R e t = 


p l U t .t D t.t 


100m 


where 


and 


p i 



p i = 




(C27) 


(C28) 
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The preceding quantities and the methods of appendix B were used to estimate tur- 
bine efficiency. As indicated, this procedure was iterated until component efficiency was 
within ±0. 0005. The remaining cycle temperatures and the cycle thermal efficiency were 
then calculated from 


(Ah), 


\ 

T 5 = T 4 + — 

(C30) 

c p 


3 * T 2 - E r< T 2 - T 5> 

(C31) 

6 = T 5 + E r< T 2 - T 5> 

(C32) 


and 
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The liquid radiator temperatures were determined from the assigned values of E , 


C , and 
s’ 


T 4, l “ T 3 ” 


T 3- T 4 

E„ 


(C34) 


= T 4,! +C s< T 3- T 4) 
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where 
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The required prime radiator area A was approximated from the following relation 

* j P 

from reference 11 
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(m V ; 

4a eT^ 
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< T 4,i - T s>< T 3, ; +T s> 
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where the radiator wall temperature drop has been neglected. 

The output from this procedure included all cycle temperatures, compressor and 
turbine pressures, overall geometry of the turbine and compressor, with characteristic 
numbers plus cycle thermal efficiency and prime liquid radiator area requirements. 
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APPENDIX D 


HEAT SOURCE ANALYSIS 

by John L. Klann, Leonard Soffer, and Gerald J. Barna 

Gross shielding approximations and factors affecting fuel block geometry and weight 
are presented. 


Shielding 

Approximate shield thicknesses which were calculated in the analysis are shown in 
figure 19 as a function of the initial plutonium 238 thermal inventory and two radiation 
dose rates. These thicknesses were obtained by an unpublished approximate technique 
which was intended for use in demonstrating trends. The fuel block was treated as a 
point source, and sell- shielding was neglected. The radiation dose distribution at the 
manned compartment was assumed to be 90 percent gammas and 10 percent neutrons in 
order to approximate minimum shield weight. 



Initial thermal power, kW 


Figure 19. - Approximate shield thickness as function of 
initial plutonium 238 thermal inventory and dose rate. 
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Lithium hydride was used for shielding the gamma and neutron radiation. Over the 
range of initial radioisotope thermal inventory of interest for this analysis (approximately 
10 to 40 kW), the shield thicknesses varied from about 11 to 14. 5 inches (28 to 37 cm) 
for a dose rate of 5 millirems per hour at a distance of 10 feet (3 m) and from about 17 
to 20. 5 inches (43 to 52 cm) for a dose rate of 1 millirem per hour. Only the lower dose 
rate required tungsten shielding (approximately 0. 1 in. , or 0.25 cm) for the allowed 
gamma dose. 

Figure 20 presents the side and front (corresponding to the side of the large face of 
the fuel block facing the power system) shield weights, calculated from the shield thick- 
nesses of figure 19, as functions of the two dose rates and the initial radioisotope inven- 
tory. The shield geometry was that idealized in figure 2. The fuel block was assumed 
to be a point source located 7 inches (17. 8 cm) behind the front shield. A fuel block 
specific area density (Ap/Q^ n Q ) of 0. 162 square foot per kilowatt (0. 015 m^/kW) was 
assumed. Hence, for a given thermal inventory, the fuel block frontal area was deter- 
mined. With a 1-inch (2. 54-cm) clearance between the fuel block and side shields and 
6 inches (15. 2 cm) of shield overhang beneath the fuel block, the dimensions of the four 
shield pieces and their weights were then determined. Radiation scattering was not con- 
sidered in determining the shield weight. A comparison was made between the weights 
shown in figure 20 and those calculated by more rigorous methods. The curves resulted 
in a total shield weight which was 10 to 20 percent heavier than those which were more 
exact. 



Figure 20. - Shield weights for assumed geometry as 
function of initial power and dose rate. 
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Over the range of parameters investigated in this analysis, fuel block specific area 
densities varied from 0. 22 to 0. 38 square foot per kilowatt (0. 020 to 0. 035 m^/kW). 
Hence, the shield weights were obtained indirectly from figure 20. The front shield 
weight was scaled directly as the ratio of the calculated fuel block specific area density 

p p 

to the reference value of the figure (0. 162 ft / kW, or 0. 015 mykW) and summed with 
the uncorrected combined side shield weight. 


Fuel Block 

The concept and methods of reference 2 were employed in a digital computer pro- 
gram to calculate the geometry and weight of the fuel block. A schematic diagram of the 
6 -kilowatt -electric powerplant fuel block is shown in figure 21. Nine rows of fuel 
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Exploded cross-section 

Figure 21. - Schematic diagram of fuel-block geometry. 
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capsules spanned the width of this block. Each capsule had an inner diameter of 2. 40 
inches (6. 10 cm) in which the plutonium fuel would be placed. A tantalum alloy (T-222) 
was used as the strength-bearing member. The T-222 wall thickness (0. 445 in. , or 
1. 13 cm, in this case) was determined to withstand the maximum pressure (12 500 psi, 
or 86. 1 MN/m^) at 2460° R (1367° K) from the helium generation during the isotope 
decay. A 0. 075-inch- (0. 19-cm-) thick layer of platinum -rhodium was employed to 
resist oxidation and Earth-substance corrosion. For the fuel block material, a high- 
melting-point material (over 3460° R, or 1900° K) with reasonable strength in the 
2460° R (1367° K) range was needed. This led to consideration of the refractory metals. 
Other desirable characteristics were relatively low density and high thermal conductivity. 
Niobium was selected as a compromise and because it does not form a volatile oxide. 
Although oxygen embrittlement will take place in the Earth's atmosphere, the block 
should remain intact. All radiant surfaces (niobium) of the fuel block were assumed to 
be coated with iron titanate for high emittance. 

Additional studies of the radioisotope heat source concept are reported in refer- 
ence 8. Our philosophy was that the radioisotope fuel must be contained during all cred- 
ible events. In reference 8, the block concept has been replaced by an array of individual 
capsules supported within a bank of tubes mounted on a reentry body. Resulting heat 
transfer is then by radiation from the capsules to the tube holders to the heat-source 
heat exchanger. The individual fuel capsules are then sized so that burial on reentry 
will not result in a capsule meltdown and fuel release. 

The present procedure assumed that all fuel -block radiant surfaces had an emittance 
of 0. 86 and an overall coefficient of thermal conductivity for the composite block of 
40 Btu per hour per foot per °F or 250 kilojoules per hour per meter per °K. Require- 
ments based on capsule burial were not included in determining capsule and block size. 

The weight and size trends are assumed to be representative of the general heat 
source philosophy, although the fuel block concept is not currently being pursued. 

Required surface area for thermal radiation interchange between the heat source and 
heat exchanger for the fuel block and the tube arrays of reference 8 were approximately 
the same. However, the heat source weights in reference 8 are somewhat lighter (ap- 
proximately 20 percent), because of the difference in construction. 

Unless noted otherwise, the fuel block reference conditions for this analysis are as 
follows: unaugmented (0 = 1 , no fins), single-sided thermal interchange between the 
block and heat exchanger; one layer of fuel capsules; capsule void volumes sized for an 

O 

internal pressure of 10 000 psi (69 MN/m ) at approximately 130 years (ref. 2). 

Figure 22 presents the variations of fuel-block surface area and weight over a range 
of turbine inlet temperatures for the reference conditions and explores some possible 
alternate limitations. For each case in this figure, the fuel block reached the maximum 
temperature of 2460° R (1367° K). Three cases are shown with limiting "normal" heat 
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Figure 22. - Effect of turbine inlet temperature on fuel-block 
surface area and weight. Gross shaft power, 6.875 kilo- 
watts; capsule internal pressure, 10 000 psi (69 MN/m^). 
(All other initial reference conditions of table I apply. ) 
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transfer (i. e. , during radiant thermal transfer to the heat-source heat exchanger). 

These are for (1) the reference conditions (8 = 1), (2) single-sided heat transfer augment- 
ed by intermeshing fins on both the fuel -block and heat -exchanger surfaces (8 = 4), and 
(3) heat transfer from both large faces of the fuel block ("double-sided" heat transfer). 

A constant temperature drop of 50° R (28° K) was assumed between the working gas and 
the heat-source heat -exchanger external surface temperature. For normal heat-transfer 
limiting, the temperature limit was reached near the fuel -block surface facing the reentry 
body (fig. 2) for single-sided heat transfer and at the fuel -block center for double-sided 
heat transfer. Two cases are shown with limiting "emergency" heat transfer (i. e. , 
when the encapsulated energy must be capable of being rejected to the environment from 
an exposed surface). These are for (1) direct dumping from the fuel-block surface as 
would occur if the entire reentry body and fuel block were pivoted out of the spacecraft 
(0° radiation dump temperature drop), and (2) indirect dumping from an equal surface 
area but at a lower temperature, as might occur if the energy could be rejected through 
the reentry body after opening the insulated spacecraft hatch (fig. 2). 

The maximum surface area and weight were required of the fuel block for the assum- 
ed reference unaugmented conditions (solid curve, 8 = 1). As turbine inlet temperature 
is increased, cycle efficiency is increased and consequently, the thermal inventory 
required is decreased. However, the temperatures to which the fuel block radiates its 
energy also is increased and thus eventually requires an increase in surface area in 
order to maintain the maximum fuel -block temperature. The required surface area is 
relatively constant up to a turbine inlet temperature of 1900° R (1056° K) and then in- 
creases while the block weight continuously decreases over the temperature range from 
1800° to 2200° R (1000° to 1222° K). The decrease in fuel -block weight was caused by 
the decreased fuel inventory and a thinner block, in spite of the increase in surface area. 

Another limiting mode for the fuel block is the requirement of emergency heat dump 
directly from the fuel -block surface (dashed line, 0° R, radiation dump temperature 
drop). This case represents the size of fuel block that can satisfy the emergency 
requirement, and if the normal requirements could be altered, it represents the mini- 
mum size. Because the normal temperature limitation was ignored, the surface area 
and weight of the fuel block continued to decrease with increasing turbine inlet tempera- 
ture as a result of the decreasing radioisotope inventory. The two heat-transfer modes 
form a limited band of fuel -block weight. The weights shown, however, are of course 
dependent on the fuel -block concept which represents a conservative approach and per- 
haps a pessimistic weight. 

Alternate arrangements for normal heat transfer might utilize either finning to 
enhance the radiation interchange, or radiation from both large faces of the fuel block, 
or both. The solid curve (8 = 4) shows the effect of finning both fuel block and heat- 
source heat exchanger so that the effective, single -sided, heat-transfer area is increased 
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by a factor of 4. The number and size of fins for this amplification has not been analy- 
zed, nor has the mechanical complexity of the physical arrangement been determined. 
However, the assumption ( 9 = 4) reduced the single-sided, unfinned surface area require- 
ments by about 27 percent. The effect of double-sided radiation heat transfer also 
reduced the fuel -block surface-area requirements for the normal -mode limitation. The 
unaugmented double-sided heat transfer reduced the surface area by about 36 percent 
from the reference case. For the emergency -mode limitation, the effect of a 100° R 
(or 56° K) temperature drop (such as might occur from the assumed geometry of a fuel 
block placed in a reentry body) before radiation to the environment, was also deter- 
mined. This additional temperature drop required the fuel -block surface area to 
increase about 7 percent over the unrestricted emergency -mode case. 




Figure 23. - Effect of cycle thermal efficiency on fuel-block 
surface area and weight. Gross shaft power, 6.875 kilo- 
watts; turbine inlet temperature, 1960 °R (1089° K); 
capsule internal pressure, 10 000 psi (69 MN/m 2 ). 
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Since each mode in figure 22 ignored the temperature limitation of the other mode, 
a composite variation of fuel -block surface area requirements with turbine inlet temper- 
ature might be required, depending on the chosen alternatives, to avoid exceeding 
2460° R (1367° K) in both these modes. For example, if the ground rules were double- 
sided, unaugmented, normal heat transfer and emergency heat transfer with a 100° R 
(56° K) temperature drop before radiation dumping, the variation of surface area with 
turbine inlet temperature would follow the dashed 100° R (56° K) line up to a temper- 
ature of 1990° R (1106° K), where a minimum surface area of 4. 2 square feet (0. 39 m^) 
would occur, and after which would follow the solid double-sided heat-transfer line. 

Although these alternatives yielded smaller surface areas and weights (which would 
also reflect in less shield weight), they were not pursued, and the single -sided, un- 
augmented, normal heat-transfer-limitedblocks were used in the text. There is an 
advantage in utilizing a large block limited for normal heat transfer; that is, during the 
emergency mode, which may be the greater portion of its life if it is not recovered, the 
maximum fuel block temperature is considerably less than 2460° R (1367° K). This, of 
course, would favor most material requirements for long-term containment. 

Figure 23 presents the variations of fuel-block surface area and weight with cycle 
efficiency for the same limiting modes in the preceding figure. For each limiting mode, 
the variation of surface area with cycle efficiency was nearly linear. Hence, the specific 
area density for a fixed turbine inlet temperature and the block reference condition were 
constant. In all cases, increasing cycle efficiency decreased fuel -block surface area 
and weight. For the reference case (6 = 1, solid line), both surface area and weight 
decreased by about 28 percent from a cycle efficiency of 0. 24 to 0. 34. 

Typical variations of fuel-block size and weight with design maximum internal pres- 
sure are shown in figure 24. Increasing the design maximum internal fuel capsule pres- 
sure continuously decreases the fuel capsule inner diameter for a fixed inventory. As the 
capsules become smaller, the overall fuel -block thickness decreases until the increased 
wall thickness required by the higher pressure overrides the decrease in capsule inner 
diameter. In figure 24, this minimum block thickness occurred at a maximum pressure 

p 

of 10 000 psi (69 MN/m ). Required fuel -block surface area continuously decreased over 

p 

the range in design pressure from 4000 to 20 000 psi (27 to 140 MN/m ), while fuel -block 

p 

weight exhibited a minimum of 7500 psi (52 MN/m^). Since the fuel -block surface area 

p 

was still decreasing above 7500 psi (52 MN/m ), the combined weight of fuel block and 

shield may still decrease above this pressure. 

The relative weights of fuel block, shield, and their sum are shown in figure 25 as 

a function of the design maximum internal fuel capsule pressure. Because of the smaller 

fuel -block face area for a fixed inventory, the shield weight continuously decreased with 

increasing design pressure. Hence, the opposing effects of fuel -block weight and shield 

2 

weight exhibited an insensitive minimum at about 12 000 psi (86 MN/m ). 
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Capsule internal pressure, psi 


3 6 9 12xl0 7 

Capsule internal pressure, N/m 2 

Figure 24. - Typical fuel-block variations with internal 
fuel capsule pressure. Gross shaft power, 6.875 
kilowatts; turbine inlet temperature, 2060° R 
(1144° K); cycle efficiency, 0.314; fin area- 
multiplication factor, 1. 
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Internal fuel capsule pressure, N/rrr 

Figure 25. - Effect of design internal fuel capsule pressure on relative heat 
source weights. Gross shaft power, 6.875 kilowatts; dose rate at 10 feet 
(3 m), 5 millirems per hour. 




Figure 26. - Typical effects of finning on fuel-block size. 
Gross shaft power, 6.875 kilowatts; turbine inlet tem- 
perature, 2060° R (1144° K); cycle efficiency, 0.314. 
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Effects of varying the number of layers of fuel capsules were analyzed for a fixed 
thermal inventory. Increasing the number of layers with the capsules in a hexagonal 
array decreased the fuel capsule inner diameter but increased the block thickness. The 
increase in layers increased the number of fuel capsules, while the required surface 
area was essentially unchanged. Fuel -block weight, however, increased by about 12 per- 
cent from one layer of capsules to two. 

Typical effects of finning on fuel -block size are shown in figure 26 for a continuous 
variation of the fin area -multiplication factor from 1 to almost 6. Required surface area 
was decreased by about 30 percent over this range, however, about 60 percent of the 
reduction occurred going from no fins to an area factor of 2. Hence, a relatively small 
amount of finning appears to go a long way in reducing the area requirements. However, 
both fuel -block thickness and capsule inner diameter required increases due to finning. 
The finned concept seems to warrant further detailed investigation. 
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APPENDIX E 


HEAT-EXCHANGER ANALYSIS 

by John L. Klann, Gabriel N. Kaykaty, Paul T. Kerwin, and Dari D. Bien 

Brief descriptions of the calculations and typical results are presented for each of 
the cycle heat-transfer components. For each of these components, the calculation 
employed digital computer programs and assumed particular core geometries. Further- 
more, a section is included on the allocation of component pressure drops. 


Heat-Source Heat Exchanger 

This component was assumed to be a simple, once-through, plate-fin heat exchanger 
and was matched to the radiation area of the fuel block and the cycle requirements. The 
large surface of the heat exchanger was assumed to be square. For a given set of 
requirements, the product of the number of fins per foot and the fin thickness was varied 
over an appropriate range. At each value of this product, the fin height was determined 
by satisfying the heat load and an assigned pressure drop. Enough solutions were obtained 
to indicate both the minimum volume and minimum weight heat exchangers. 

Figure 27 presents the heat-source heat -exchanger weights and volumes as a func- 
tion of percent pressure drop Ap/p and system power level. The solid lines show the 
constant power minimum volume heat exchangers which were used for system perfor- 
mance calculations. The dashed lines indicate some typical minimum weight results at 
constant power. As system power level was decreased, heat-exchanger heat load and 
pressure level also were decreased. These effects are opposing; that is, for the same 
pressure level, decreasing heat load would decrease weight, while for the same heat 
load, decreasing pressure level would increase weight. The net effect shown in fig- 
ure 27 was the relatively small change in weight and volumes among the three power 
levels. Above about a 3-percent pressure drop, all three power levels exhibited the 
same rates of change of slightly decreasing weights and volumes with increased pressure 
drop. At lower values of pressure drop, the 2-kilowatt heat -exchanger weights became 
heavier than either the 6- or 10-kilowatt solutions. 

For comparison, the minimum -weight solution at the 2 -kilowatt power level and the 
volume for the 10-kilowatt-power -level minimum -weight solution are presented in 
figure 27. With increased pressure drop, the difference between minimum weight and vol- 
ume decreases. In the range of pressure drop of interest (2 to 3 percent), for 2 kilowatts 
system power, the minimum -weight solution was about 30 pounds (13. 6 kg) lighter than 
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Figure 27. - Effect of pressure drop on reference system heat-source heat 
exchangers. 

the minimum -volume solution, while at the 10-kilowatt level, the minimum -volume 

O 

solution was about 0. 6 cubic feet (0. 017 m ) smaller than the minimum -weight solution. 
Since this exchanger was included in the shielded volume, the relatively small increase 
in heat-exchanger weight probably would be justified by use of the thinnest heat exchanger. 

The relative effect of working fluid and molecular weight on the heat-source heat- 
exchanger weight (minimum -weight solutions) was explored and is presented in figure 28. 
A helium -xenon mixture at a molecular weight of 83. 8 is the reference, and other, work- 
ing fluids are compared with it. Differences among the working fluids were smallest at 
the 10-percent pressure drop where the use of krypton would increase the heat -exchanger 
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Figure 28. - Effect of working fluid on relative weight of 
heat-source heat exchanger (minimum-weight 
solutions). 


weight 68 percent. However, in the range of interest (e.g. , at a 3-percent pressure 
drop), changing from a helium -xenon mixture, at the respective molecular weight, to 
either argon or krypton would about double the heat-exchanger weight. A similar rela- 
tive weight relation was noted for the recuperator. 


Recuperator 

The recuperator in this analysis was assumed to have the same construction as the 
recuperator in reference 1. A rectangular, offset fin, counterflow core with triangular 
end sections was used. The high-pressure-side and low-pressure-side core geometries 
were those of figures 10-58 and 10-59, respectively, of reference 12 with the splitter 
plates removed. The core face area was assumed to be square. For a given set of cycle 
conditions, the length for a desired heat load (including' effects of axial conduction) and 
effectiveness was calculated as a function of the face area. Pressure drops were then 
determined from manifold to manifold for each side including contraction, expansion, 
turning, and friction losses. Total recuperator percent pressure drops were used as the 
sum of the high- and low-pressure-side percent pressure drops. 

Figure 29 presents the variations of total recuperator weight and core volume with 
total percent pressure drop as a function of assigned effectiveness. These are typical 
results and were those used in the construction of figures 13 and 14. A dashed line indi- 
cates the values of pressure drop which were assigned to the recuperator for minimum 
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(a) Recuperator weight. (b) Recuperator volume. 

Figure 29. - Effect of pressure drop on recuperator weight and core volume. Gross shaft power, 6.875 kilowatts. (All other 
initial reference conditions of table I apply. ) 


combined heat -exchanger weight. The assigned pressure drops increased from 3. 8 per- 
cent at an effectiveness of 0. 875 to 6. 9 percent at an effectiveness of 0. 99. Both weights 
and volumes became large as full recuperation was approached. For values of effective- 
ness between 0. 875 and 0. 95, the changes in weight variations with pressure drop were 
nearly the same. However, at higher values of effectiveness, the variations in weight 
were a function of the level of recuperation. 
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Heat-Sink Heat Exchanger 


1 


The heat-sink heat exchanger was assumed to be a multipass, cross-counterflow 
device with a finned tube core. The core was that of figure 10-84 in reference 12. The 
exchanger was sized to fit directly on the appropriate end section of the recuperator. 
Then, the n um ber of passes was increased from 1 until the desired effectiveness (0. 950 
±0. 001) was attained. Associated core pressure drops were calculated, and as an ap- 



Combined percent pressure drop 

Figure 30. - Variations of reference system recuperator and heat-source 
heat-exchanger with combined pressure drop. 
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proximate allowance for manifolding, the gas-side percent pressure drop was doubled. 
Since the mating of the recuperator and heat-sink heat exchanger was assumed, the 
weights and percent pressure drops of these components were combined for the analysis. 

Figure 30 presents the heat-sink heat-exchanger weights, the combined recuperator 
and heat-sink heat- exchanger weights, and the heat-sink heat -exchanger percent pres- 
sure drop as functions of the combined pressure drop and system power level. At 
2 kilowatts, the heat-sink heat exchanger was at most 45 percent of the combined weight, 
whereas at 6 and 10 kilowatts the heat-sink heat exchanger was less than 23 percent of 
the combined weight. The reduction of power level results in the conflicting effects of 
reduced heat load and pressure level. The 6-kilowatt combined heat exchanger above a 
combined pressure drop of 2 percent exhibited the smallest weights. 


Pressure Drop Distributions 

The criterion employed to determine the gas pressure drop split among the heat 
exchangers was that which resulted in minimum combined weight of the recuperator, 
heat-sink heat exchanger, and minimum -volume heat-source heat exchanger. This 
criterion is satisfied when the slopes of the weight against pressure drop curves for each 
of the three components are equal. 

Figure 31 presents the pressure drop distributions used in this analysis. Percent 
pressure drops are shown for the recuperator, heat-source and heat-sink heat 
exchangers along with the system total for the investigations of the effects of recuperator 
effectiveness, loss pressure ratio, and system power level. The piping percent pres- 
sure drop was assumed to be 1 percent for the investigations of recuperator effectiveness 
and loss pressure ratio, while 1. 5 percent was assumed for the final reference cases. 

For all cases, the recuperator needed the largest pressure drop, while the heat-sink heat 
exchanger needed the smallest. As recuperator effectiveness was increased, the pres- 
sure drop allocated to it also increased with the necessity of decreasing allocations to 
the remaining components. For increasing loss pressure ratio, all distributions 
decreased linearly. For a total system loss pressure ratio at the 0.92 level from 2 to 
10 kilowatts, the recuperator and heat-sink heat exchanger percent pressure drops 
decreased, while the heat-source heat exchanger percent pressure drop increased. 


Radiator 

A cylindrical, unoriented radiator with longitudinal tubes and circumferential headers 
on a 260-inch- (6.6-meter-) diameter vehicle was assumed. Two tubes (one redundant) 
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were assumed for each fin, and calculations were made for a probability of success of 
0. 995 against meteoroid penetration for a 5-year mission. FC-75 (a fluorochemical) was 
assumed to be typical of the class of coolants suitable for this application, and its heat- 
transfer properties were used in the calculations. The computer program of refer- 
ence 18, modified for cylindrical geometry, was employed. Whipple's 1963A meteoroid 
flux model was assumed along with a meteoroid density of 0. 2 gram per cubic centimeter 
and an average velocity of 65 000 feet per second (19 800 m/sec). Heat rejection was 
assumed to occur only from the outside cylindrical surface, and the inside surface armor 
was reduced because of the bumper effect by the method of reference 19. Liners inside 
the tubing were assumed with a resulting spall factor of 1.4. 

Typical results for the radiators are shown in figure 32. The effects of system 
power level on radiator area and weight are shown. The curves shown are the loci of 
minimum radiator weights as a function of surface area; that is, the number of tubes on 
the fixed diameter was varied until the combined effectiveness of tubes and fins resulted 
in a minimum weight. For each power level, the minimum value of required radiator 
area appears to approach a value of 65 square feet per kilowatt (6. 04 m /kW) of system 



Figure 32. - Effect of power level on reference system radiator area and weight. 
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output. Minimum radiator weights, on the other hand, were reached at a value of about 

O 

125 square feet per kilowatt electric output (11.6 m /kWe). 

In those parts of the analysis where fixed radiator areas were assumed, radiator 
weight variations are the result of reoptimization of the number of tubes and fin thickness 
for the changes in radiator conditions with the parameter in question. 
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APPENDIX F 


SYSTEM EFFICIENCY CALCULATIONS 


Overall system efficiency used in this analysis was defined as 


p 

a, net 


(FI) 


where 


Vr L " (Q i. ,I i ) 


P a, net P a " P k P p 


(F2) 

(F3) 


The raw alternator terminal power was defined as 


P s r? P 
a 'a, em sh, net 


(F4) 


where 


P. , = P , _ L -L - L. - L- 
sh, net sh w s j t 


(F5) 


where 


P sh 7? cy^in 


(F6) 


Hence, 


73 (P . - L - L_ - L. - L,) - P, - P 
'a, enr sh w s j t; k ; 

sy 1. 048 (Q in + h t ) 


(F7) 


Allowances were made for heat-source thermal loss L^, system control power P k> and 
radiator pump power P while estimates were made for alternator electromagnetic effi- 
ciency 7 ? a em , windage loss L w , shaft seal loss L g , and bearing friction losses Lj 
and L t . ’ 
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Gross shaft power P gh was assumed to be directly proportional to cycle mass flow 
rate. Corrections of turbine and compressor efficiency due to variations in Reynolds 
number (eq. (B7)) were small and hence were neglected. The required cycle thermal 
input Q in for a fixed gross shaft power was determined from the calculated cycle effi- 
ciency (eq. (F6)). The heat-source insulation loss at the end of the 5-year mission 
was assumed to be 5 percent of the cycle thermal input. The gross initial heat-source 
energy Qj n5 g then followed from equation (F2) where the factor 1. 048 allowed for plu- 
tonium 238 decay over the 5-year mission and an assumed 90-day launch delay. 
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Figure 33. - Effect of alternator type, power level, and cavity pressure on efficiency and 
windage losses. 
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Results of an unpublished alternator sizing study which were employed in the calcu- 
lations of this analysis are shown in figure 33. Alternator windage loss as a function of 
cavity pressure and electromagnetic efficiency as a function of gross alternator power 
output are presented for two rotational speeds and two alternator types. Alternator and 
bearing cavity pressures were assumed to be 75 percent of the compressor discharge 
pressure. For rotational speeds of 24 000 and 36 000 rpm (2510 and 3770 rad/sec), the 
rotor diameters for either type of alternator were conservatively estimated as 4. 2 and 
3. 6 inches (10.7 and 9. 14 cm), respectively. Although the Lundell-type alternator had 
lower electromagnetic efficiencies than the homopolar-inductor type, its windage losses 
were considerably smaller than those of the homopolar-inductor type. On the basis of 
overall performance, the Lundell-type alternator was preferred and the 36 000-rpm 
(3770-rad/sec) characteristics were used in the analysis. 

The alternator windage loss L w from figure 33 was used as one of the terms in 
equation (F5) to correct the gross shaft power to the net value available to the alternator. 
Other shaft windage losses were assumed to be included in the turbine and compressor 
efficiencies. The seal leakage loss L g assumed that 1 percent of the compressor flow 
was lost to the turbine to perform useful work. 

Journal bearing analysis was based on the three-tilting-pad theory of reference 20, 
while the thrust bearing analysis employed the spiral -groove theory of reference 21. 

The journal bearing loads were assumed to consist of the rotor weight and an 
arbitrary constant load of 3. 3 pounds or 14. 7 newtons (an attempt to account for other 
unbalanced forces). The rotor weight was scaled from the results of reference 3 by 
assuming a fixed weight of 8. 3 pounds (37 N) and varying the remaining weight of 
6. 2 pounds (27. 6 N) by the square of the compressor tip diameter divided by 3. 0 inches 
or 7.6 centimeters (the compressor tip diameter of ref. 3). The total journal bearing 
load was then assumed to be equally split between two sets of bearings. The shaft radius 
was determined from the 2-kilowatt system bearing cavity pressure so that either a 
maximum bearing unit loading of 4 psi (27. 6 kN/m ) or a maximum total bearing load 
coefficient of 1. 0 was attained. Then, the bearing design curves of reference 20 for a 
shoe arc length of 94. 5°, a ratio of shaft radius to shoe length of 0. 606, and a dimension- 
less pivot location of 0. 667 were employed. A preload factor of 0. 6 and a built-in shoe 
clearance to shaft radius ratio of 0. 001 were assumed. The padf riction factor was nearly 
constant over a wide range of shoe eccentricity ratio; hence, the frictional forces (cal- 
culated by eq. 57 of ref. 20) and losses were assumed constant and independent of both 
load and pressure level. 

The calculated thrust bearing losses assumed that, at a gross shaft power of 6. 875 
kilowatts, a net thrust load of 10 pounds (44. 5 N) could be attained. This arbitrary load 
was varied directly with gross shaft power. The inner shaft radius was assumed to be 
the same as the shaft radius at the journal bearing locations. For the thrust bearing, an 
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inner- to outer-radius ratio 0. 5 and a ratio of bearing clearance to outer radius of 0. 001 
was assumed. A. double-acting bearing was used with a total bearing gap of 0. 007 inch 
(0. 0178 cm). Equations (Dl), (D2), and (D4) of reference 21 were used to calculate both 
the loaded and unloaded bearing face clearances and friction losses. 

The final system penalties were allowances for control power requirements and 
radiator pump power. Control power was assumed to be 4 percent of the alternator out- 
put, while pump power was assumed to be about 10 watts per kilowatt of thermal energy 
to be rejected. Subtraction of these two losses from the gross alternator output yielded 
the net alternator power. 
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